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The effects of synoptic waves on the dynamics of planetary waves are 
investigated using normal mode analysis. Initialized analyses of the Navy Operational 
Global Atmospheric Prediction System (NOGAPS) for 19 days between January and 
April 1986 are projected onto the normal modes of a linearized version of the model. 
For each analysis, the different terms (adiabatic nonlinear, linear and diabatic) which 
affect the time tendency of planetary-scale modes are determined by a one-time step 
integration of the NOGAPS model. The effect of synoptic scales on planetary scales is 
determined by computing the difference between the adiabatic nonlinear term 
computed from the NOGAPS analyses and analyses for the same period that have 
been spectrally filtered to remove most of the svnoptic-scale waves. The energy 
tendency due to the nonlinear adiabatic term and the synoptic-scale contribution to 
this term are also computed. It 1s shown that the synoptic-scale contribution to the 
adiabatic nonlinear term and the time tendency of planetarv-scale modes can be a very 
large percentage of these terms. 

By eliminating momentum advections in the model and computing the adiabatic 
nonlinear term for the filtered and unfiltered analvses, the relative importance of 
interactions through mass field interactions or momentum field interactions are 
determined. It is shown that svnoptic-scale interactions which affect the planetarv- 
scale barotropic modes are primarily through the momentum advections, while mass 
and momentum interactions are possible for the baroclinic modes. The importance of 
mass field interaction generally increases as the vertical scale of the wave decreases. 

Because of the importance of svnoptic waves to the dynamics and energetics of 
planetary waves, errors in the forecasts of planetarv waves may in part be due to the 


svnoptic-scale forecast errors. 
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I. INTRODUCTION 


Atmospheric predictability studies (Lorenz, 1969) indicate that each scale of 
motion has different predictability, e.g. small-scale motions are theoretically predictable 
to an hour, synoptic scales to a few days and planetary scales to a few weeks. 
However, a number of studies (Lambert and Merilees, 1978; Baumhefner and Downey, 
1978; and Morse, 1983) have shown that synoptic-scale motions, rather than planetarv- 
scale motions, are the most accurately forecast by numerical models. Although 
planetary wave errors have been reduced in current operational models (Wallace et al., 
1983; Bettge and Baumhefner, 1984), planetary waves are still not forecast as accurately 
as cyclone scale waves. This result does not seem to be dependent on the type of 
model used. Baumhefner and Downey examined a number of different models and 
found the same error pattern. The error in forecasting planetary-scale waves might not 
seem to be that important given that the error in the synoptic scales is smaller and that 
the synoptic scales of motion produce most of the dav-to-day weather changes. The 
importance of accurate planetary-scale predictions is clear when one considers medium- 
to long-range forecasts (up to 10 days). Since planetary waves often steer the smaller 
synoptic disturbances, an improved planetary planetary-scale forecast would 
presumably lead to an improvement on the synoptic scale. 

There are at least two reasons why planetary waves are not forecast as well as 
theory suggests. One reason is that the dynamics for planetary waves are poorly 
understood. Planetary waves consist of a quasi-stationary component and a generally 
smaller transient component. The quasi-stationary component is thought to be a 
response to forcing by a combination of topography, differential heating due to land- 
sea differences and possibly scale interactions. It would seem reasonable to suspect 
that errors in the quasi-stationary component might be due to inadequacies in the 
model forcing. Another source of error for planetary waves is spurious excitation of 
transient planetary waves (Daley et al., 1981; Somerville, 1980; Lambert and Merilees, 
1978; Roads and Somerville, 1982). These spurious waves are thought to be the result 
of errors in initial data or an inadequate model domain. Such waves are a major 
contributor to planetary wave error in the first 24 to 48 hours of a forecast (Daley et 


al., 1981). These transient waves have often been examined in terms of the normal 


Ml 


modes of the linearized set of equations describing oscillations of a stratified resting 
atmosphere on a spherical earth. The mode that has been most often identified as 
being spuriously excited is the so-called external Rossby mode. This mode has a large 
phase speed and a period of approximately 5 days. This fact has led some to refer to 
this mode as the 5-day wave. Daley et al. (1981) showed that large-scale external 
Rossby modes are excited when a hemispheric domain is used rather than a global 
domain, and when bad or poorly analyzed tropical data are used. Daley also showed 
that application of the nonlinear normal mode initialization procedure of Machenhauer 
(1977) to the large-scale external Rossby modes improved the forecast in a root mean 
squared sense for the cases examined. However, he properly notes that there is no 
evidence to show that the Machenhauer balance condition, which seems to be 
appropriate for the gravity modes, is applicable to the fast Rossby modes. 

Very little is known about these Rossby modes in the atmosphere or in numerical 
prediction models. The best observational evidence for the existence of these normal 
modes is from the studies of Madden and Julian (1972) and Ahlquist (1982). Applying 
time series analysis techniques to station pressure data, Madden and Julian were able 
to identify westward propagating 5-day waves. The observed characteristics of these 
waves were shown to be not inconsistent with those of a planetary-scale Rossbv wave. 
Ahlquist projected 1200 consecutive days of twice daily National Meteorological 
Center tropospheric analvses of velocity and geopotential onto three-dimensional, 
normal mode Rossby wave structures. Through spectral analyses of these time series, 
Ahlquist was able to identify 14 planetary-scale, normal mode waves. By contrast, a 
rather large amount of research (Dickenson and Williamson, 1972; Williamson, 1976; 
NIachenhauer, 1977; Errico, 1984) has been done to determine the nature of the 
dynamic balance of the gravity modes. It is from these investigations that the 
nonlinear normal mode initialization procedure was developed. The dynamic balance 
of the gravity modes has been studied primarily by long-term model integrations. 

A third and as yet unexamined reason for errors in numerical predictions of 
planetary waves is that these errors are due to errors in the smaller (cyclone) scales. 
The dynamics of planetary waves may be such that nonlinear interactions from smaller 
scales are important. If this is the case, then errors in cyclone scales would lead to 
errors in the planetary scales. The planetary scale prediction errors in numerical 
models cannot be fully understood or corrected until the importance of nonlinear scale 


interactions are determined. 
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It is the general hypothesis of this study that nonlinear interactions with cyclone 
waves are an important factor in the dynamics of planetary waves. The purpose of this 
study is twofold: 


1. Determine the importance of these nonlinear interactions using the normal 
mode analysis technique of Errico (1984); and 


2. Identify some of the mechanisms of these nonlinear interactions. 

A number of studies have examined nonlinear interactions and their importance 
to the maintenance of planetary waves. One of the earliest studies to illustrate how 
nonlinear interactions could maintain large-scale, quasi-permanent flow was by 
Saltzman (1959). Using a highly simplified barotropic model with an idealized flow 
that crudely simulated a Northern Hemisphere winter 500 mb flow pattern, Saltzman 
found that there was a substantial energy transfer from the synoptic waves to the 

planetary waves. Saltzman and Fleisher (1960) used 500 mb data to show that in the 
mean there was a net kinetic energy transfer from synoptic scales to planetary scales. 
A more recent study bv Kao and Lee (1977) showed that the primary contribution of 
nonlinear interactions to the energy transfer is essentially through the interactions of 
the slowly moving waves, the stationary long waves and zonal mean flow. Saltzman 
(1970) gives a review of the major studies that have used Fourier analysis to identify 
nonlinear interactions. Gall et al. (1979) used a simplified general circulation model to 
demonstrate that the initial development of the ultralong waves from a zonal mean 
basic state can be forced by the interaction between the cyclone waves and the basic 
flow. Gall argued that the ultralong waves were forced mainly by planetary scale 
Variations in the meridional heat flux convergence of the higher wavenumber modes, 
which produces a positive correlation between planetary scale upward motion and 
temperature. The principal kinetic energy source for the planetary waves was the 
conversion of wave available potential energy to wave Kinetic energy at a given 
wavenumber. A more recent study by Young and Villere (1985) confirmed in part 
Gall’s results, but also showed that direct transfer of kinetic energy from intermediate 
scales to planetary scale was of equal importance. In both of these studies, the 
nonlinear transfer of potential energy was not computed directly. This transfer was 
implied by showing that the conversion of eddy available potential energy to eddy 
kinetic energy was much greater when nonlinear interactions were allowed. 

The analysis of the scale interactions in this study are done using Errico’s normal 


mode analysis procedure. This analysis procedure requires the use of a numerical 
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model and its associated normal modes. The model that will be used in this study is a 
version of the Navv Operational Global Atmospheric Prediction System (NOGAPS). 
The normal modes used in the analysis procedure are derived from this model. This 
analysis method is used because it has a number of advantages over previously used 
methods such as Fourier or spherical harmonics analysis. The most important 
advantage of using normal mode analysis is-that the total effect of other scales 
(vertical, meridional and zonal) on a given mode can be determined. Another 
important advantage is that the nature of wave motions can be identified physically 
and high-frequency noise can be separated consistently from low-frequency 
meteorologically significant waves in both mass and velocity fields. Normal mode 
analysis also allows one to separate the data into different vertical, zonal and 
meridional scales while maintaining the physical nature of the data. The other 
techniques that were mentioned above do not have this advantage. While it is true 
that a combination of the Fourier and spherical harmonic could be used to decompose 
the data into different vertical, meridional and zonal scales, this representation would 
be artificial and it has the disadvantage of having no explicit relationship between 
spectral modes of mass and velocity. Finally, since divergence is significant in the 
motion of ultralong waves, it may be more appropriate to represent data as solutions 
of linearized primitive equations rather than as solutions of the non-divergent vorticity 
equation (i.e. spherical harmonics). 

There have only been a few studies in recent years (Kasahara and Puri, 1981; Ko, 
1985) that have used normal mode analysis to examine the spectral distribution of 
atmospheric energy. These studies were confined to exanuning just the spectral 
distribution of atmospheric energy and not the energy conversions between modes of 
different vertical and horizontal scales . There has been even less work in determining 
the energy exchange between the different vertical and horizontal modes. The little 
work that has been done in this area has focused on zonal mean-eddy kinetic energy 
exchange (Tanaka et al., 1986). Tanaka used normal mode analysis, to examine 25 
davs of dailv First GARP (Global Atmospherics Research Program) Global 
Experiment (FGGE) IIIb analyses from the Goddard Laboratory for Atmospheres 
(GLA) and the Geophysical Fluid Dynamics Laboratory (GFDL). Tanaka examined 
the distribution of Kinetic energv as well as the kinetic energy interactions between the 
barotropic mode and baroclinic modes of different zonal wavenumbers. The focus of 


his study was on zonal mean-eddy kinetic energy exchange. This study is unique in 


Id 


that the focus of the normal mode analysis is on the interaction between synoptic 
scales and planetary scales and that these interactions are examined in term of the total 
energy exchange (potential and kinetic) between scales. 

A description of the NOGAPS model used in this study is given in Chapter II. 
The normal modes of this model are derived in Chapter III. To assist in the 
interpretation of results from the NOGAPS model, a nonlinear scale analysis and the 
results from a simple analytic model are presented in Chapter IV. The results of the 
experiments conducted to determine the importance and mechanisms of nonlinear 
interactions are given in Chapter V. Chapter VI contains a brief summary and some 


conclusions of this work. 
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Il. THE MODEL 


The model used in this study to investigate the dynamics of planetary scale waves 
is version 2.2: of the Navy Operational Global Atmospheric Prediction System 
(NOGAPS). This model is a modified version of the general circulation model 
developed at the University of California at Los Angeles (Arakawa and Lamb, 1977). 
It is a global finite difference model and uses the primitive equations in sigma 
coordinates. The resolution of the mass variables (surface pressure and temperature) 1s 
2.4° latitude by 3.0° longitude with nine levels from 50 mb to the surface. The model 
uses scheme C staggering (Arakawa and Lamb, 1977) in the horizontal, and finite 
differencing in the vertical according to Arakawa and Suarez (1983). The horizontal 
finite differencing 1s energy conserving and it conserves enstrophy when the motion 1s 
nondivergent. The vertical differencing conserves the global mass integral of the 
potential temperature under adiabatic processes and it emplovs a local form of the 
hydrostatic equation. The time differencing is a combination of five leapfrog steps for 
each Matsuno backward step, while the heating 1s computed during a single forward 
step that precedes the Matsuno step. However, in this study only the forward time 
step of the Matsuno step is used 1n the analysis procedure. 

Phvsical processes include radiation, moist and drv convective adjustment and a 
cumulus parameterization technique (Arakawa and Schubert, 1974), which interacts 
with a bulk parameter boundary laver (Randall, 1976; Lord, 1978). 

The objective analyses of wind and geopotential are done with 4a three- 
dimensional successive corrections method that is a form of the scheme used by Barnes 
(1964). The analyses of wind and mass are done independently. Since the NOGAPS 
initialization method 1s fully described by Barker (1982). onlv a brief description of this 
method will be given here. The results of the independent wind and mass analvses are 
combined via a calculus of variation method in which the balance equation 1s used as a 


constraint. The functional 


F(4,¥) = | (6 — 3)? +B(V — V)? +2M1V- SY ‘ 
+2J(u,v) — V*¢|dA 
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is minimized and the associated Euler-Lagrange equations are solved over the entire 
earth. Here 9 is geopotential, y is the stream function, A is the horizontal area over 
which the integral is applied, 4 is a Lagrange multiplier and the symboi ~ denotes 
analyzed values. The quantity Bp can be made a function of latitude if desired to force 
more adjustment toward either one of the analysis fields in certain areas. Thus, B 
would have relatively greater weight in low latitudes where the wind is a more reliable 
parameter for analysis and prediction than the geopotential field. The above procedure 
minimizes the change to the analyzed geopotential and non-divergent wind while 
constraining these variables toward the balance condition. The problem of generating 
an appropriate divergence to go with the nondivergent winds produced by this 
balancing procedure is solved by using the forecast first-guess divergence. The problem 
of vertically inconsistent corrections is minimized by vertically coupling the variables. 
The variables are coupled before they are initialized by projecting them onto empirical 
orthogonal functions. The smoothness of the four empirical orthogonal functions used 
insures that the inconsistent vertical variations of wind or geopotential that could be 
generated by the initialization procedure are eliminated. 

The initialization procedure used in this version of the NOGAPS model initializes 
objectively analyzed correction fields for V and > rather than the updated fields. This 
procedure has the advantage of not affecting areas without new data. The resulting 
initialized corrections are interpolated to the model sigma coordinates surfaces and 
added to the first-guess forecast. This method minimizes vertical interpolation error 


and preserves the model generated first guess divergence. 
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Il. NORMAL NIODES 


Since the data used in this studv are analyzed using the normal modes of tne 
linearized equations from the NOGAPS model, it 1s important to have a fundamental 
understanding of these modes. This understanding may be gained bv examining the 
derivation and structure of the modes. The normal modes for an earlier (lower 
resolution) version of the NOGAPS model were derived bv Barker (1982). He followed 
the method of Temperton (1977) and Temperton and Willamson (1981). Three 
different sets of modes are derived for the version of the NOGAPS model used in this 
study. Thev are obtained below bv separation of variables. Except where noted, their 
derivation closely follows that of Barker (1982). Each set of modes is examined for 
consistency with the model as well as for simularities to modes derived by other 
authors. [t will be shown that modes derived with a model top at 0 mb are best suited 


for use in this study. 


A. VERTICAL MODES 


The linearized governing equations used by Temperton and Williamson (1981) 





are 
oY fk xV + V(RT inp, +4) =Q, 3.1) 
= Heo (se 
SALOME 9 hy, Sat LOY (3:3 
at 
eer acetal (3.4) 


Here the vertical discretization is taken into account by writing the equation in vector 
form. Thus, V is the vector form of the wind, T is the perturbation temperature, T is 
the rest-state temperature, p, is surface pressure, $ is the perturbation geopotential, 
V°V is divergence and 9, is terrain geopotential. t and G are linearized matrix 


Operators and I! is a vector. Q,, Qz and Q, are the nonlinear components of their 


{8 


respective equations. All of the above operators and vectors are defined in Temperton 
and Williamson (1981). Their exact form is not important. However, it 1s important 
to note that the entries that make up these matrices are constant and that o is defined 


as p/P. 
Barker approximated (3.1) - (3.4) using the o system of the NOGAPS model as 








+ fkx V+ V(oralnt +4) =Q, (3.3) 
Bae (VV) =Q, 3 (3.6) 
ot 
oS" +I" (V-V) =Q, a 
ot 
d¢=¢,+GT (3.8) 
where 
a eet (3.9) 
T 
and 
T=), — Dy. (3.10) 


p, is the pressure at the top of the model atmosphere and @ is specific volume. t and 
G are Once again matrix operators that are similar to those in (3.2)-(3.4). Thev are 
consistent with the Arakawa and Suarez (1983) vertical finite difference scheme used in 
the NOGAPS model and are given in Appendix A. A different treatment of (3.4) will 
be discussed later. 

Following Temperton and Williamson (1981), Barker (1982) defined a vector h 


Whose horizontal derivative represents the pressure gradient force, 


gh=¢+o607aln7 (3.11) 


This definition allows the determination of a single equation for mass by operating on 





(3.6) with G, multiplying (3.7) bva7a@, and then adding the resulting two equations to 
obtain 


oh 
In; +C(V-V) =Q, (eb2) 


where 


C =Gr+orall’ (3.13) 


and 


n = GQ7 + Q,RT. _—. (3.14) 


Although equation set (3.12) is vertically coupled, it can be transformed by 
separation of variables into a set that is not coupled. This is done through the 


diagonalization of the matrix C by: 


a | ~ a 
ECE = gD, (3.00 
where the matrix E contains the eigenvectors of the matrix C as columns and the 


diagonal matrix D contains the eigenvalues of C. Defining the vertical transforms of h 


and V as 
h=E-th (3.17) 
produces the uncoupled equation set 
av 3 ee 
a tik xV +oVb=Q, (3.18) 
dh aN 
95, + D(V-V) =Qa, (3.19) 


Where Q,. and Q, are the transforms of Q,, and Q,, respectively. These equations are 
uncoupled except through the nonlinear terms on the right sides of the equations. The 
independent variables in (3.18) and (3.19) are the coefficients of the vertical modes (the 
eigenvectors contained in E). There are as many modes as there are tevels in the 
model. 

The vertical modes derived using the linearized NOGAPS equations given by 
(3.5)-(3.8) for T equal to (218, 218, 218, 228, 241, 254, 267, 276, 283)°X are show meu 
Fig. 3.1. For comparison, the verucal modes derived by Temperton and Willamson 
(1981) for T equal to (229, 209, 218, 237. 256, 268, 277, 283, 285)°K are"shown mies 
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Figure 3.1. Vertical modes for the NOGAPS model with nine levels, T equal to (218 
218, 218, 228, 241, 254, 267, 276, 283)°K anda top at 50 mb. 


b 


al 


how 


Verret? @ tt 


— 





Figure 3.2. Vertical modes from the model used by Temperton and Williamson. The 
model top is at 0 mb and T is equal to (229, 209, 218, 237, 256, 268, 277, 283, 233)4 
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Table 1. Equivalent depths (m) of the vertical modes of the NOGAPS and 
Temperton and Williamson (T & W) models. 
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mode index 
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3.2. The eigenvalues (equivalent depths) of the vertical modes for this linearization of 
the NOGAPS model and those derived by Temperton and Williamson (1981) are given 
in Table 1. These NOGAPS modes are qualitatively very similar to the modes derived 
by Temperton and Williamson (1981). The general structure (except the sign of the 
4th-6th modes are opposite from Temperton and Williamson, but this an arbitrary 
choice) and the levels at which the crossings occur are similar. However, there are 
some noticeable differences between the two sets of modes, with the most significant 
differences occuring in the gravest modes. The NOGAPS barotropic mode (mode one) 
lacks the characteristic peak at the top of the atmosphere and in fact decreases near 
the top. The first two baroclinic NOGAPS modes have zero crossings at levels lower 
than for the modes derived by Temperton and Williamson (1981). Most of these 
differences can be attributed to the differences between the NOGAPS model and the 
model used by Temperton and Williamson (1981). One difference between models is 
the location of the model levels. One of the major differences. in terms of vertical 
Structure, 1s the location of the model top. The model top in the NOGAPS model is at 
50 mb while the top in the model used by Temperton and Williamson (1981) is at O 
mb. Barker (1982) also found that the equivalent depths are sensitive to the location 
of the model top. Changing the model top from 50 mb to U mb increased the 
equivalent depth of the external mode from 7874 m to 966V0 m. Barker found that a 
consequence of specifving the model top at 50 mb ts that all the equivalent depths are 
smaller than if the top was at VU mb. 

Although most of the difference between the two sets of modes can be attributed 


to the difference in the model tops, part of this difference is due to an inconsistent 


lineanzation of the hydrostatic equation by Barker. The discretized, unlinearized 
hydrostatic equation for the NOGAPS model ts 


@¢=¢%,+ GT, (3.20) 
where the matrix operator G is similar to the one given in Appendix A, except that the 
full values and not the mean state values are used in the entnes of G. The G matrix 
defined by Temperton and Williamson (1981) is constant. However, the G matnx 
defined above is not constant. The choice of a model top at a nonzero pressure level 
makes G a function of surface pressure (m in this case), which was not taken into 
account by Barker (1982) in his linearization of the hydrostatic equation (3.8). 

A second set of NOGAPS vertical modes that includes the etfect of u variable G 


matrix will now be derived. A more consistent linearization of (3.20) vields 


¢ = $,+GT'+G'T, GoM 
where 
G = 4G] CG. (3.229 


This linearization of (3.20) has an effect on the matnx C {rom which the vertical modes 


are derived. Consider the momentum equation 


OV 

By tik x V tonaVine + Vo =Qy. (3.23) 
Inserting (3.20) in (3.23) 

OV 

By + fk x Vt+onaVinn + Vd, + VGT = Qy. (3.24) 


The matrix G ts a function of m. which is a function of the horizontal coordinates. 
Expanding VGT one obtains 

OV 

a + fk x V toraVina + Vg, +(VG)T+GVT=Q,. 3.25) 


Following the suggestion of Rosmond (1986) VG can be written with the chain rule as 


dG 
Wa NS ln (3.26) 
Using (3.26), (3.25) can be linearized as 
av' _dG— a 
Tag dis 2S NY + (ora +7 T)V Inn + V$, + GVT = Ge eyra 
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The derivative of G with respect to m can be determined analvtically since G can be 
Written as a function of x. The linearization of the remaining governing equations is 


unchanged. To obtain a single equation for mass one must define h as 


ace a 
gh = (ora + 7—_—T) ln op, Gr - (3.28) 
Using the above definition of h, the matnx C would be defined as 
= _dG 
= Gr + (oma + 7 


—— T)MI*. (3.29) 


The matrix C defined in (3.29) is clearly different from the C defined in (3.13). If 
G were not a function of surface pressure (the model top was at zero), then (3.28) and 
(3.29) would reduce to (3.11) and (3.13) respectively. Since the matrix C in (3.29) 1s 
different from that defined in (3.13), the eigenvectors will be different. However, the 
overall similarities between the vertical eigenvectors and eigenvalues of the NOGAPS 
model (as derived bv Barker, 1982) and those of the model used bv Temperton and 


Williamson (1981) indicate that the effect of the additional term 
dG — 

tat (3.30) 
in (3.29) is small. The vertical modes obtained from the matrix C in (3.29) are shown 
in Fig. 3.3 and the equivalent depths are given in Table 2. The term (3.30) has been 
approximated bv a centered finite difference as 

_ G(# + At) + G(w — AZ) Fr 

ino. oar ei 
Where Az equals | mb. The major difference between these modes and those derived 
previously occurs in the barotropic mode. The characteristic peak that is absent in 
mode one in Fig. 3.1 is now present. Also, the equivalent depth associated with this 
mode has increased from 7829 m to 8101 m. The more consistent linearization of the 
hydrostatic equation leads to a set of vertical modes that are more similar to those 
derived by Temperton and Williamson (1981). 

Although the vertical modes derived with a variable G matrix are more simular 
(than modes derived without this effect) to modes derived by other authors, they do 
have some disadvantages. They require the computation and storage of an additional 
term. In addition, the placement of a model top at 0 mb is clearly a more realistic 
condition than having the model top at a finite pressure level. The closer the basic 


State is to the real atmosphere, the better the linear approximation becomes. Thus, it 


Table 2. Equivalent depths (m) of the vertical modes of the NOGAPS model that 
includes the effect of a variable G matrix. 
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would seem that better results would be obtained at less computational and storage 
expense if the vertical modes were derived with a model top at 0 mb. 

Consequently, the modes used in this study are derived from a basic state in 
which the top of the atmosphere is assumed to be at 0 mb. The equivalent depths for 
these modes are given in Table 3 and their structure is given in Fig. 3.4. The modes 
derived with the model top at 0 mb have equivalent depths that are much larger than 
the NOGAPS modes derived using Barker’s linearization or those derived with a 
variable G matrix. The 0 mb top modes have Zero crossing at higher levels than the 
two other sets of NOGAPS modes and their structure is more similar to the modes 


derived bv Temperton and Willamson. 


B. HORIZONTAL MODES 

The details of the determination of the horizontal modes of the NOGAPS model 
have been given by Barker (1982). Onlv the general method for determining the modes 
will be given here. 

The solutions of the uncoupled equation set (3.18 - 3.19 ) for each equivalent 
depth Dy gives the horizontal modes of the model. Equations (3.18) and (3.19) can be 


Written in finite difference form as 
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Fig 3.3. As in Figure 3.1, except modes are derived including the effect of a variabie G 


Matrix. 


Table 3. Equivalent depths (m) of the vertical modes of the NOGAPS model with 
model top at 0 mb. 


Vertical Equivalent depth 
mode index 
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= (65%)i, , (aM)5) 5 (3.34) 
ehs + DUG tp abd! Phe : 


A vertical mode index (2) is assumed for D and each variable. The finite difference 


operators are 





(6.7). = ae = hd a (3.35) 

——z Te 4 +7ZI,_1 

(T), = - (3.36) 
Vy Saar 

(7). = (oe: (3.500) 


The other variables are defined as follows: a is the earth's radius, AA is the longitudinal 
grid interval, A® is the latitudinal grid interval, u and v are the east and north 
components of the wind, respectively, 1 and } are the longitudinal and latitudinal 
indexes, respectively and fp is cos 9. 

Special definitions of the Coriolis term are used to keep the matrix operator of 
(3.32)-(3.34) symmetric. A symmetric matrix insures that the corresponding 
eigenvectors are orthogonal, which allows determination of the inverse of the 
eigenvector matrix by simply taking the transpose. To achieve symmetry, the Coniolis 


term in (3.32) is replacedsby: 
za =< as —- 


2p, (3.38) 
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Figure 3.4. Same as Figure 3.1, except model top is at 0 mb. 
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and the Coriolis term in (3.33) is replaced by 
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; (3.39) 
where 
ee aebieek a5 ses 
: cos(4*) (3.40) 
and 7 
ie Spee a afi 
cos(“") (3.41) 


These definitions correspond to a potential enstrophy conserving finite difference 
scheme as derived by Temperton and Williamson (1979). 
A dynamical state vector is defined as 
TAC 3) 
MOnses es S| eae) (3.4m) 
A(A;, 65, 2) 


where y represents grid point values of the vertical mode coefficients. These values can 


be expanded into Fourier modes by 
(oar 
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The inverse of (3.43) 1s 
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ailleagete te) =]! AACR we ee 
7 6) (3.44) 
Assuming a wave solution of the form in (3.43) allows (3.32) - (3.34) to be written as 
Ou; ., (7th) oe | gik’h, _s 
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vy _ sin(44) 
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mes (3.49) 





r(m) = cos( 


AA EX 
J = —) — tsin(— 
(m) = cos(—") ~ ésin(>). “5 
After applving the symmetric operator 
wan) 0. 0 
vue i a fens 
0 0 (g/D.)? 
tory such that 
go = Se4 (3.52) 
allows Eqs. (3.45) - (3.47) to be written in matrix form as: 
: Oy’ ~) “C\ ET! 
—ae + Ly’ = -:QH’. (3.53) 


The matrix Q is diagonal and positive definite with entries related to the cos 0). is a 
matrix depending on the actual finite differences used 1n the model, and H is a matrix 
representing the latitudinal variation of the corresponding nonlinear terms transformed 

in the same wav as u, v and h are to give y 
To complete the computation of the horizontal modes, the matrix equation (3.53) 


1s rescaled using 
a— Q: +, H =Q7H’, L=Q°?LQ"?. (3.54) 


If Y contains the eigenvectors (normal modes) of L, then (3.53) can be written using 
(3.54) as 


.d 2 : 
Creve Gale) + Y-*LY(Y~?4) =-rY-'H. 
The identity 
A= Y ‘LY. (3.56) 


where A is a diagonal matrix containing the eigenvalues of L, makes it possible to 
rewrite (3.55) as 


dC 
ay =-7AC+Yr, (3.57) 


where 


a 


C=ay- "4. (3.58) 
The nonlinear term is now r and the mode frequencies are A. 
The components of C, denoted here by C(m,n,£), are referred to as the coefficient 
of the normal mode. Equation (3.57) can be written in component from as 
dC (m,n, £) 
Bt (3.59) 
The elements of C are functions of the vertical, zonal and meridional mode numbers @, 


= -twC(m,n,£) + r(m,n, £). 


m and n, respectively. The mode number n is a measure of the number of zeros 
between the poles. For each m, 2, and n there are three equations for C: One for an 
eastward propagating gravity wave (EG); one for a westward propagating gravity 
(WG); wave and one for a westward propagating Rossby wave (R). The symmetric 
modes (u and h are symmetric about the equator and v is antisymmetric) consist of the 
odd indexed (n=1,3,5..) Rossby. modes and even indexed gravity modes. The 
anuusvmmetric modes (u and h antisymmetric about the equator and v 1s symmetric) 
consist of the even indexed Rossby modes and the odd indexed gravitv modes. These 
coefficients are the amplitudes of the various modes required to represent a particular 
atmospheric state. Corresponding to each horizonal mode ts a natural frequency (@) 
that is determined as an eigenvalue of the system. The frequencies of various modes 
corresponding to those given by Dickenson and Wilhamson (1972), Temperton and 
Williamson (1981) and Barker (1982) are given in Tables 4 and 5. The frequencies of 
all the modes are very similar, with most of small differences being due to the different 
horizontal resolutions used in each model. The structures for a few selected middes are 
given in Figs. 3.5 and 3.6. The structures of these modes are also very similar to those 
derived by other authors (Temperton and Williamson, 1981: Kasahara,1976; Dickenson 
and Wilhamson, 1972). Anvy results obtained here using the modes of the NOGAPS 


model should be similar to results obtained using modes of other models. 


C. EXPANSION OF DATA INTO NORMAL MODES 

The amplitude of a given mode is determined bv expanding grid point values into 
normal modes. The first step in the expansion process is to remove the mean state, 
then combine the thermodynamic variables into one variable, the equivalent 
geopotential (gh). The data are expanded into vertical modes using (3.16) and (3.17), 
then into Fourier modes using (3.44) and lastly scaled according to (3.51). The 
symmetric and antisymmetric components are found bv averaging or diflerencing the 
values from the two hemispheres. Finally, each scaled Fourier mode of each vertical 


mode is expanded into meridional modes bv (3.58). 
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Table 4. Frequencies (s7!) of the Rossby modes for the models used by Temperton 
and Williamson (1981) (T&W), Dickenson and Williamson (1972) (D&W), for the 
NOGAPS model with model top at 50 mb (B) and for the NOGAPS model with 
model top at 0 mb (M). D=10 km for (T&W) and (D&W), D=8,101 m for (B) 
and D= 9,682 m for (M). Horizontal grid intervals are specified in degrees. 


n T&W, 10 D&W, 2.5 Bt2. 433 M2.4 x 3 
0 6.11x107> 6.14x107> 5.93x107> 6.14x107° 
1 1.44x1075 1.45x107> 1.32x107> 1.43x107° 
2 8.46x107© 8.73x107& 98.16x107© 8.68x107® 
3. 5.72x107© 5.87x107© 5.61x107© 5.86x107° 
4 3.98x107© 4.17x107© 4.05x107© 4.18x107° 
5 2.87x107& 3.08x107© 3.03x107© 3.11x107° 
6 2.14x107& 2.36x107© 2.34x107° 2.39x107® 
7 1.63x107& 1.86x107© 1.86x107© 1.89x107° 
8 1.27x107& 1.49x107& 1.51x107& 1.52x107° 
9 1.01x107& 1.22x107© 1.24x107& 1.25x107§ = | 
10 8.10x1077 1.02x107& 1.04x107& 1.05x107® 
11 6.62x1077 8.58x1077 8.83x1077 8.89x107’ 
12 5.52x1077 7.30x1077 7.58x1077? 7.62x107’ 
13. 4.70x107’ 6.27x107? 6.57x107’ 6.60x107/ 
14 4.11x1077 5.431077 5.74x1077 5.77x1077 
15 3.75x1077 4.73x1077 5.06x1077? 5.08x107’ 
16 3.13x107’ 4.14x1077 4.49x1077 4.50x1077 
| 


For each 2 and zonal wavenumber m .2 0, there are 111 modes in the NOGAPS 
model for the symmetric case and 112 for the antisymmetric case. These modes mav 
be divided into three sets (sets of 37 each for the symmetric case, two sets of 37 and 
one set of 38 for the antisymmetric case). The modes associated with the 37 largest 
negative eigenvalues are usuallv referred to as westward-gravitational (WG) modes. 
The 37 modes whose eigenvalues are positive are refered to as eastward-gravitational 
(EG) modes. The remaining modes, all with negative or zero values are referred to as 


Rossby or rotational (R). In a linearized model, these modes describe westward- and 


Table 5. As in Table 4, except for eastward-gravitational modes. 


n T&W, 10 D&W, 2.5 Bw -X oO M2.4 x3 
0 -5.44x107>2 <-5.38x107> -4.71x1072 -5.29x107> 
1 -1.31x107* -1.30x107% -1.18x1074 -1.28x107* 
2 -1.87x107* -1.86x107* -1.69x107* -1.84x1074 
3. -2.35x107* -2.36x107* -2.13x107* -2.33x1074 
4 -2.79x107* -2.83x107* -2.55x107* -2.80x1074 
5 ©3.22x107* -3.29x107* <2.97x107" =3. 27m 
6 -3.63x107% -3.75x107* -3.38x107% -3.73x1074 
7 #4.01x107* %<4.21x107* -3.80x10"=" ==4. Zone 
8 -4.36x107* -4.66x107* -4.22x107* -1.51x1074 
9 -4.69x107* -5.10x107* -4.63x107* -4.66x107* 
10 -4.98x107* -5.54x107* -5.05x107* -5.13x107* 
11 -5.23x107* -5.97x107* -5.47x107* -5.60x107% 
12. -5.44x107* -6.38x107* -5.88x107* -6.06x107* 
13. -5.61x107* -6.79x107% -6.29x107% -6.52x1074 
14 -5.72x107* -7.19x107* -6.70x107* -6.98x1074 
15 <-5.94x107* -7.57x107* -7.11x107= "=7. 44x10 
16 -5.94x107* -7.94x107* -7.51x107* -7.89x107- 
eastward-propagating gravitv waves, and westward-propagating Rossby waves, 


respectivelv. Consistent with the above categorization. Kelvin waves are EG modes 
and mixed Rossby-gravitvy waves are R modes. For each 2, there are 113 zonaliv 
symmetric (m=() modes in the NOGAPS model. The modes may be divided among 
the tvpes R. WG and EG. All of the R modes for m=V are stationary in the linearized 
model. 

In the real atmosphere or in nonlinear models, the normal modes from the 
linearized equations are no longer independent solutions and their behavior is not 
necessarily wave-like. However. each mode may be considered as describing a 


particular dynamical structure. The behavior of these modes is altered by the inclusion 
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Figure 3.5. Structure of the u component (dotted), v component (dashed), and the h 
component for selected rotational modes (R) of the model in this study. I is the 
Zonal wavenumber and N is the meridional mode number. 
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Fig 3.6. As Figure 3.5, except for selected gravitational modes (EG, WG). 


of those terms/physical processes that are absent in the linearized models: nonlinear 
advective terms and diabatic terms (which include friction). These terms couple the 
behavior of the modes in a nonlinear model. The eigenvalue associated with each 
mode in a nonlinear model is not the only frequency the mode may have, but it may be 
considered as the mode’s natural or resonant frequency. Whether a mode with that 
frequency is observed depends on the strength of the nonlinear and diabatic forcing. 
For the modes of equivalent depth De the structures of the modes have been 
normalized so that the sum of the kinetic plus available potential energies per unit 
mass of the corresponding fields is given by 
1 * 

E.= 5 >» OER an 
Where C,, is the amplitude of the mode designated by the index n and the asterisk 
denotes a complex conjugate. The sum is over all modes of a given equivalent depth 
Dp (denoted by the set Sp). For each value of £, each mode contributes independently 
to Ep because the modes are mutually orthogonal for each ee 

The vertical modes of the NOGAPS model are not orthogonal. Therefore, it is 


not true that 


9 
ra \ 
dX : (3.61) 


is the total kinetic plus available potential energy of all the modes. However, E is 
positive definite and increases in value as the Ep increase. Although its is not the 


precise energy, it does have utility. 


D. NORMAL MODE ANALYSIS 

In the nonlinear, diabatic, discrete-tume NOGAPS model the prognostic equation 
of a coefficient mav be written as 

a =-wt, +N, +Q, (3.62) 

where A C,,'At is the discrete time tendency of a parucular vertical. zonal. meridional 
EG. WG, or R mode, -10C, represents the linear terms in the model, i=V1 , Sa 
represents the adiabatic nonlinear terms, and Q,, 18 a sum of diabatic terms. The focus 
of this study is on the prognostic equation for the R modes and the nonlinear 
interactions which may occur through the adiabatic nonlinear term. 

The procedure for determining the different terms in (3.62) for the R modes is 


relativity simple. The linear term for a particular mode n is determined by multipling 


oF 


the coefficient of a mode by its natural frequency. The coefficient for a mode is 
determined from initialized NOGAPS analyses using the procedure described in section 
C of this chapter. The discrete time tendency is determined from a one-time step 
integration of the NOGAPS model. First the coefficients for the modes are determined 
at the initial time (t). Then the model is integrated over one time step (At) and the 
coefficients att + At are determined. The discrete time tendency is then given by 
C,,(t + At) — C, (t) 

At (3.63) 
where At is equal to four nunutes. The determination of the time tendency was found 
not to be sensitive to the time step used. In addition, for the Rossbv modes, one 
would not expect AC/At to vary greatly from one time step to another. This in fact 
has been demonstrated bv Errico (1984). However, for the gravitational modes, 
AC/At mav vary largely in the initial stages of a model integration. This vanation of 
AC/At for the gravitation mode will depend on the initialization method used. For the 
Rossby modes, the procedure for determining the time tendency should be adequate for 
the purposes of this study. 

The nonlinear adiabatic term is determined bv subtracting the linear term from 
the discrete time tendency that has been determined from an adiabatic one-time step 
integration. Errico (1981) was able to derive explicit expressions for the Ns of a 
simple, nonlinear, f-plane, primitive equation, two laver model. For this simple model, 
the Ns can be written in terms of the sums of quadratic functions of the modal 
amplitudes. These sums represent interactions between rotational (geostrophic)}, 
rotational and gravitational (ageostrophic) and gravitational modes of different 
wavenumbers. The complex NOGAPS model does not lend itself to such a simple 
decomposition. Thus the nonlinear adiabatic term is determined numerically using the 
NOGAPS model and the procedure described above. 

Using a numerical procedure, Errico (1984) was able to decompose the adiabatic 
nonlinear terms of the modal equations of a gobal pmmitive equation spectral model 
into three different groups: NW(Ri. N,(G) and N,(R*G). The first depends only on 
the rotational mode coefficients. the second only on the gravitational mode coefficients, 
and the iast only on the sums of products of each tvpe of mode with the other. 
Because of the highly divergent nature of the gravational modes. Ny(R*G). may be 
likened to interactions between the divergent component of the motion and the 
rotational component. For a simple f-plane model, \,, (R*G) would represent terms 


such as the advection of momentum bv the divergent wind. 
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The adiabatic nonlinear terms for the modal equations of the NOGAPS model 
are not decomposed into N,(R), N,(G), and N(R*G). However, Errico (1984) has 
shown for the long time solutions to the NCAR Community Climate Model, N(R) is 
an order of magnitude larger than N,(R*G) and N,(G) for all but the smallest vertical 
scale rotational modes (£=8,9). For the smallest vertical scale R modes N,,(R*G) ~ 
N(R). For rotational modes, one would expect that if nonlinear interactions occur 
through the adiabatic nonlinear term that they would to a first order approximation, be 
due to interactions between rotational modes (except for the smallest scale vertical 
~ modes i.e. 2= 8,9). 

The diabatic term is determined by subtracting the discrete time tendency 
determined from the adiabatic one-time step integration from the discrete time 
tendency determined from a one-time step integration which included diabatic 
parameterizations (e.g. cumulus, sensible heating parameterizations etc). 

In addition to allowing the examination of possible nonlinear interactions 
through the adiabatic nonlinear term, (3.62) can be used to gain insight into the 
general dynamical nature of a given mode. This insight may be gained by comparing 
the relative magnitudes of the different terms. For example, it is possible to determine 
if a mode 1s balanced, 1e., the time tendency of the mode is small compared to the 
other terms. [fa mode is not balanced, it may be possible to ascertain which terms are 


most important in determining the time tendency of a given mode. 


IV. ASCALE ANALYSIS AND A SIMPLE ANALYTIC MODEL 


To interpret the results from the complex NOGAPS model, one must first have a 
good understanding of the linear and nonlinear dynamics of planetarv scale modes. To 
highhght the difference between the results presented in the next chapter and those 
predicted by linear theorv, a linear and nonlinear scale analvsis of the shallow water 
equations is presented. A simpler analytic model will also be developed in this section. 
This simple model will be used to show what terms are represented by the adiabatic 


nonlinear term and how svnoptic waves mav affect planetary waves. 


A. SCALE ANALYSIS 

Some important dynamical features of planetary scale Rossby modes can be 
ulustrated bv a simple scale analvsis of the shallow water equations with a vanable 
Coriolis parameter. The shallow water equations are adequate for this purpose 
because the primitive equations can be represented through the normal mode as a set 
of shallow water equations for each equivalent depth (assuming V is not a function of 


the vertical coordinate). The shallow water equations to be scale analyzed are 


OV = 


and 

oe +V-Vd+6V-V+iovV-V. 

ot (4.2) 
The nonlinear terms have been retained to show how nonlinear interactions mav have 
a significant affect on planetary scale motions. In a strictly shallow water context. o 
would be considered the mean height of the fluid and would be fixed. In a normal 
mode context, $ is the equivalent depth and is different for each vertical mode 
considered. Note that 9 is the departure of the geopotential from >. 

Two length scales are chosen to demonstrate how synoptic scale motions (L>) 
affect planetary scale Rossby motions (Lj). In addition, two different ume scales, an 
advective scale (IT =L,V) and one appropriate for a free planetarv scale Rossby wave, 
will be used. The latter time scale is determined by L,/C, Where C = BO: f* is the 
phase speed of the fastest Rossby mode. Using two length scales and the time scale 


associated with this fast Rossby mode makes the scale analysis presented in this section 
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different from those of other authors (Burger, 1958; Haltiner and Williams, 1980) who 
used an advective time scale and a single planetary length scale. The scales used in the 


following analyses are: 


f= 10-*s**, L,=10'm, ,i2 = 108m, V=10ms"' (4.3) 


1. Single length scale analysis 
Using the planetary length scale given above and assuming an advective time 


scale, the equation of motion (4.1), with the exception of the pressure gradient force, 1s 


scaled as 
V 
Boy VV +V6+fkxVi=0 | 
at (4.4) 
2 ie 
ma vw & 
i Ly, L, 
or after multiplication by L, 
y2 y? d fV 
so that the values are 
10> 1102 d 10* 


If the equation is to be balanced, the pressure gradient force must balance the coriolis 
force i.e. the scaling is geostrophic, so that > scales to 107, 
Using geostrophic scaling for 0, an advective time scale and a single planetary 


length scale, the geopotential equation (4.2) 1s 


oP 4 V-V6+8V-V46V-V (4.5) 
vo =e ov va 
Dy, 1, Ly, Ly 
and after muluplication by L, 
Vo Vo OV Vo 


For the above analysis, @ is the scale value of > and it Represents tie Scale of the 
equivalent depth. The appropriate value of ® for the fast time scaling is 10°m*s7~ (this 
is the approximate value for the equivalent depth of the barotropic mode), Using this 
value of ® and fast time scale in (4.5) reveals that the time tendencv and the 
divergence terms are of the same order. For the slow time scaling, if one assumes 


@=10° then the first-order equation for (4.5) is V°V=0. However, most of the 
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planetary scale modes with a natural period (i.e, the period determined from the 
natural frequency) that 1s approximately equal to the slow time scale do not have an 
equivalent depth equal to 10°. A more appropriate equivalent depth for these modes is 
oye OP. Using the slow time scale and the more appropriate value of ® shows all the 
terms in (4.5) to be the same magnitude. The results obtained using the advective time 
scaling and the more appropriate value for the equivalent depth are similar to the 
results obtained by other authors (Burger, 1958; Haltiner and Williams, 1980). 

The first-order system obtained using the fast time scale is quite different from 
that obtained using the slow time scale. The fast time scale svstem is linear, the winds 
are geostrophic and the time variation is given solely by the divergence term in (4.5). 
The slow time scale system is nonlinear, the winds are also geostrophic, but the time 
Variation 1s determined by advection of the mass field as well as by the divergence term. 

2. Multi-length scale analysis 

To accomplish a multi-length scale analysis, the synoptic length scale (L4) is 
used whenever a derivative is taken in the nonlinear terms. Planetary scaling (Lj) is 
used in the other terms. A more formal procedure to determine nonlinear interactions 
would be to transform the equations using the appropriate transform based on the 
geometry: Fourier transforms for cartesian geometry or spherical harmonics for 
spherical geometry. After transformation of the equations, the nonlinear terms would 
appear as interaction coefficients that represent explicitly the scale interactions. For 
scale analysis purposes, this interaction can be illustrated by using only two length 
scales (L,, L5) in the nonlinear terms. Using multiple length scales in the scaling of 
(4.1) yields results similar to those obtained using a single length scale except the 
advection term is an order of magnitude larger. However, this term ts still an order of 
magnitude smaller than the largest term in the equation. 

Using multiple length scales in the scaling of (4.2) gives the same results as 
using a single length scale if one assumes geostrophic synoptic scaling for 9 ie 
O=fVL,. The overall effect of including nonlinear interactions for the fast Rossby 
waves 1s small. The momentum advection term in (4.1) is an order of magnitude larger 
but it is still an order of magnitude smaller than the largest term in (4.1). Because of 
the assumption of geostrophic synoptic scaling for 9 in the nonlinear analysis of (4.5), 
there is no difference between the linear and nonlinear scaling of this equation. For 
the fast time scaling, the time tendency in (4.5) is still driven by the divergence term 


and for the slow scaling (with = 10*) all the terms are still the same magnitude. 
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It would appear from the analysis using the slow scaling and $=10* that 
nonlinear interactions from the advection of the synoptic scale momentum fields may 
not be as important as the nonlinear advection of the mass field. This result is 
consistent with geostrophic adjustment theory, which indicates that the atmospheric 
response to a mass perturbation for scales larger than the Rossby radius Y@ /f is greater 
than for a given wind perturbation. The planetary length scale used in this analysis is 
an order of magnitude larger than the Rossby radius. Since the Rossby radius 
decreases as the equivalent depth (>) decreases, one would expect mass perturbations 
to be even more important for the smaller equivalent depth planetary scale baroclinic 
modes. 

Additional information on the dynamics of planetary scale Rossby modes can 
be obtained bv scale analyzing the quasi-geostrophic potential vorticity equation. This 


equation can be derived from (4.1) and (4.2). The potential vorticity equation for this 














model is 
o¢ =f OD i 
—_—_ sa Wey pe VY +vf =0 
Mag +’ YS 3Y¥9 
Lt d rE d L, 
Cc fCcL, JopetA fV I, 1 
Lif jae ¢ 


Where the nonlinear scaling for this equation is based on fast time scaling using 


6=10°. For this scaling, these nondimensional quantities have the values: 
C 1 fcl, L,V 


l 1 —S— 


my 10 ¢ fL @ 10 


Note that the nonlinear interaction through the advection of vorticity is now part of 











the first-order svstem. This implies that nonlinear momentum interactions could be 
important. Again. this result is consistent with geostrophic adjustment theorv because 
for these modes 9=10° the Rossby radius is approximately the same order of 
magnitude as the planetary length scale. 

If the nonlinear terms are dropped and the geostrophic relation for V is 


inserted, (4.6) can be written in cartesian coordinates as 


—— + = 0, 
got faz (4.7) 


Insert the wave solution ¢ = exp <j x-ct) into (4.7) then a 
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= of? (4.8) 

This agrees with the value of C that is used in this scale analysis. However, the effect 
of the nonlinear terms may be such that this phase speed is not observed for this mode, 
so that it is probably an overestimate of the speed of this mode. One final point to be 
made from this analysis is that the importance of the mass field interactions will 


increase as the equivalent depth (>) decreases. 


B. ASIMPLER ANALYTIC MODEL 

The normal modes of the NOGAPS model are quite complex and are not easily 
interpreted. A simpler analytic model will be developed in this section to illustrate 
_ what terms, effects are represented by the nonlinear term of a mode. Using this model, 
it will be shown how waves of a particular vertical and horizontal scale can affect 
waves of different scales. The results from this quasi-geostrophic model will be used *> 
interpret the more complex NOGAPS model. 


Following Haltiner and Williams (1980), the quasi-geostrophic equations are 





OG 0 —-Z7\ _ 
at Te -V(¢+ Boy) — foe” sale Z) =O (4.9) 
8 a¢' a¢ : 

as my 2S Z)\Z =0 | 

at OZ Me ag th ) (4.10) 
V’¢' — fo =0 (4.11) 
where 

@ .@ 
ee ey (4.12) 
Vy =foikx Vo=kx Vu (4.13) 
5 = f0'V7¢ (4.14) 
006 -— H*g,g 10T ; 

i Z SS = —_—_ — — —}, a. 

(2) = Sa lagt *e) = = (cS Haz? (4.15) 
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The vertical coordinate (Z) of this system is defined and related to the more familiar 
pressure coordinate system in Appendix B. The vertical velocity of this system is given 
by Z. Now making the Boussinesq approximation ( e = constant ) and assuming 
static stability (I-) is constant, then (4.9)-(4.11) may be written 


A _, 2, , ow OZ 
ae p+kx Vy~y-VVip+ 7, bo — faz =0 - (4.16) 
Ad aw op Te 


If the above is linearized about a basic state of rest and Z = 0 at the surface (Z = 0) 


and at the top of the atmosphere (Z = Z7). then the vertical modes for this system are 


ex Z 
®,(Z) = Ay cos = 
Lr 





(4.18) 
The derivation of these modes is given in Appendix B. To transform (4.16) and (4.17) 


let 
b= D> velzyy, tet (4.19) 
a = : -g ig. 
a= Z(x,y,t)e #7 
a (4.20) 
where ‘Pp=W_p and Zp = -Zp Substituting the above expressions into (4.16) and 


(4.17), muluplying bv e(XhaZiZy) and integrating from Z=0 to Z=Z7 gives the 


following transformed equations 


O O ; 
a Ve + Sk x Voy -VV2 5 + aT + fitZ, =0 eels 
: 
Owe ) ae 
£— + €—7)k x Vo, -Vi_; + —-Z = 
at ; mdi ¥7° VP ~ act ip “sie: (4.22) 
Where the boundary conditions require Zp = 0 if 2 = 0 (£=0 in this system 


corresponds to £=1 for the NOGAPS mode, 1.., the barotropic mode). The second 
terms in (4.21) and (4.22) represent the nonlinear interactions between vertical modes J 
and 2-s that will affect a given mode 2. The interaction term in (4.21) represents 
interaction between the momentum fields, while the interaction between the mass fields 


(temperature advection) is represented by the interaction term in (4.22). 
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It will be useful to look at two cases. For the first case 2 = 0 (barotropic case), 


so that (4.21) and (4.22), reduce to one equation 

< V7 Wo +S ok x Vd, -VV?p_; + we =0 

: (4.23) 
As expected, this equation closely resembles the barotropic quasi-geostrophic vorticity 
equation, except that the simple advection term is represented as the nonlinear 
interaction term. The only interactions that are allowed come from the momentum 
advection terms. The interaction through the mass advection term is eliminated 
because of the boundary conditions Z = 0 at the top and bottom of the model. The 
boundary condition at the bottom 1s not the same as in the NOGAPS model where the 
vertical motion is not necessarily zero at the bottom. One consequence of the 
boundarv condition Z = 0 at the top and bottom of the simple model is that the 
vertical structure of the barotropic mode 1s strictly constant. Also, since the equivalent 
depth is proportional to 1,2, the equivalent depth of the barotropic mode is infinite. If 
the equivalent depth (9) is infinite then from (4.6) it can be seen that mass advection 
term 1s zero. The equivalent depth associated with the barotropic mode of the 
NOGAPS model is not infinite. Thus, one would not expect that interactions through 
the mass advection term would be completely eliminated. However, the scale analysis 
of (4.6) indicated that for an equivalent depth that is approximately equal to that of 
the barotropic mode of the NOGAPS model, that interactions through this term would 
be small. It remains to be seen if the results for the barotropic mode of this simple 
model applv to results obtained with the more complex NOGAPS model. 
The second case 1s just the general baroclinic case where £2 = 0. For example, if 

& = 3 then (4.21) and (4.22) become 





a ; 2 Avs a 
7 ws +)Ukx Vo; VV vs-3 + ZT Bo + fo3Ze = 0 (4.24) 
Pj 
Ov, . ey 
= (s-3) b; > Vibs-; Ae (4.25) 


Thus, interaction mav occur for baroclinic modes through the temperature advection 
term and through the momentum advection term. For mode number 3 to be affected 
by nonlinear interactions, mode W, will interact with Wo Ws with Wa ete. 

Horizonal wave interaction will now be examined using Fourier transforms. The 
following derivation closelv follows the method used in Lorenz (1960) and Haltiner and 


Wilhams (1980). If we assume periodic boundarv conditions in both x and y 
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(4.26) 


ve(z + =, y+ pot) = vel Y, t) 


27 
k 
then the an variations of Wp and Zp can be represented using orthogonal basis 


Z(z+ — += =, t) = Z,(z, y, t) (4.27) 


functions of the form: 
t(mkz+ngy) 


Pmn iz, y) =e (4.28) 


Using the above, y and Z are represented as 


(x,y, Z,t) = 2 2 2, Comnlt n(theimas tere! te (4.29) 
(x,y, Z,t) t) = 2 2 2, Pemn(t zl er er he * Or) : i (4,30) 


where 2 and g are — Perici in the x and v directions respectively. Therefore, (4.21) 





and (4.22) can be transformed to 








OCm.t = et t +e si eV = EN El x NEC, C rues 
Ot 4.31 
(4.31) 
+5 Doc 
ACM. (eae) be 
- = B 2 —-H)-(M-H)k-HxMCyu-u.;Cuye_; (432) 
Jj 
Tr 
ae 7PM. fy 
0 


where M = mkt + ngs, \S’=¥ and H is a dummy index. The interaction terms in 
the above equations represent (aneenauens between vertical as well as homzontal 
modes. The horizontal modes interact in the same manner as do the vertical modes. 
That 1s, zonal wavenumber 7 mav interact with zonal wavenumber 8 to affect zonal! 
wavenumber |. 

Eliminating Day p between (4.31) and (4.32) gives the spectral form of the 


potential vorticity equation 


OC'm.e tM KkBo Cry e ieee A, 
hs BPS? ly ae oe —_ 
J 
where 


47 


Ai =(M-H)-(M-H)k-HxM 
and 


A, = (M-M) + 


The spectral form of the quasi-geostrophic vertical motion equation is obtained by 


fe 
. 


eliminating OCyy p/dt between (4.31) and (4.32) to give 








Due = —— Ta emy tea) — ale a + Lom. ngCn.e- 5 (4.34) 
where 
PRE ih 
A 
>= (Sea tz 


Using (4.33) and (4.34) some analogies to the normal modes of the NOGAPS can now 
be made. 
The spectrai form of the quasi-geostrophic potenual vorticity equation (4.53) 1s 


of the form’ ds (5.62) 1e 





NC. 
ra — —1wC,, aE Na (4.35) 
where 
oa (4.36) 
and 
foe Ay 
Nn, = (1+ -——) — Cnn Cue 5 
r a Be Az (4.37) 


Thus (4.33) mav be likened to the prognostic equation for the rotational modes of the 
more complex NOGAPS model. For this simple model the adiabatic nonlinear term 
(N,) 1s written as the sum of quadratic products of the modal coefficients and it 
represents interactions between the rotational modes of different wavenumbers. Note 
that there is no prognosuc equauon for Dyy p. Day p may be likened to the coefficient 
of a gravitational mode for the more complex NOGAPS model. The equation for 
Dy 2 is purely diagnostic and the amplitude of Dny e is determined solelv bv the 
rotauonal components of the motion (this is consistent with the initial assumption of 


quasi-geostrophic conditions). A final point to be made from (4.33) and (4.34) is that 
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N,, for the quasi-geostrophic equations does not contain interactions between the 
divergent and rotational components of the motion. 
It is desirable to form the kinetic and potential energy equations for these modes. 


The kinetic energy for this system is given by 


rd 2a/ ir l 
k= [ /, af Vod- Ve dZdzrdy = 5 oe tM MOm eC- me,—¢-(4.38) 
oO 0 0 | ¢ M —_ 


ao 


using the above and (4.16), the kinetic energy equation for this system may be written 
as 
OK _1 —H)-(M—H)k-H x M@y Cue; 
SFE Comet LM B) (Me Mans CHa, 
rd M 3 
(4.39) 
+ folCm,-eDm.,e 
= ~HCy,-e-3 
Ot \_(-M - 8) .(—M — H)k H x ( M)C_m-uHCu,-¢-3 
j HH 


— folC_m,--D-m-e- 


A potential energy equation may be formed by assuming the potential energy is 


proportional to 
| lt ) ) fol o, \ Oe ere 
T 9 9 | | 
to oM (4.40) 


Using the above and (4.22), the potential energy equation for this system may be 








Written as oP _ 
at 
2 
3 , = C_m-c), >) (€-s)(M-H-(M- B)k-H x MCm- uj Cue-s 
cS aa (4.41) 


+ folCmu,-2Dm.e 
2 


(Ke C | 
C | ie —-M —H.(—-M-—H)k-H x (-M)C_m-u Teen 
r mt De 2 7)( 





- 


— folC_m,-eD-m-e 
Based on (4.39) and (4.41), the amount of kinetic or potential energy transferred to a 
given mode Cnr 2 is proportional to the amplitude of that mode times the pair of 


interacting modes CuL-H jCH 2-5: Also note that the potential energy as defined in 
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The scale analysis and simple analytic model results presented in this chapter 
indicate how the dynamics of planetary-scale waves could be affected by synoptic-scale 
waves and how this effect might vary with vertical mode. However, the simple system 
of equations used here do not address the other important aspects of the dynamics of 
planetary-scale waves. For example, the importance of forcing by diabatic processes 
and nonlinear interactions which involve the divergent component of the wind are not 
addressed by the quasi-geostrophic equations. For any given time, these neglected 
effects could be small, but in a climatic sense (time mean) they may be important. The 
adiabatic nonlinear term (N_,) for the rotational modes of the NOGAPS model will 
include interactions between different scales of the highly divergent gravitational modes 
and the rotational modes. However, for the long-time solution to the NCAR 
Community Climate Model, Errico (1984) has shown that the adiabatic nonlinear term 
that represents interactions between the gravitational modes and rotational modes is 
generally an order of magnitude smaller than the adiabatic term that contains only 
interactions between rotational modes (except for the shallowest vertical modes 
£=8,9). 

While the results of the simple analytic model do not contain all the possible 
interactions between modes, one would expect the results to be valid as a first order 
approximation to the more complex NOGAPS model. In addition, the possible 
importance of diabatic process and forcing by the divergent part of the motion do not 
negate the possible importance of synoptic-scale interactions through the adiabatic 


nonlinear term. 
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V. RESULTS 


The first experiment of this study is designed to determine the dynamical 
importance of synoptic-scale nonlinear interactions on planetary scale Rossby modes. 
In the first part of this experiment, 19 days of 12UTC initialized NOGAPS analyses 
taken every five days from 16 January 1986 to 16 April 1986 are projected onto the 
normal modes of the model and the magnitudes of the terms given in (3.62) are 
determined. In addition, the energy spectrum as a function of meridional mode 
number for each vertical mode is determined. Mean magnitudes of the terms for zonal 
wavenumbers 1, 2 and 3, are obtained by averaging the real and imaginarv parts of the 
respective terms over the period and then taking the magnitude: of that average. An 
alternative way to compute these would be to take the magnitude and then average. 
The latter method has been used by most other authors (Errico, 1984; Kasahara and 
Puri, 1981). Similarities between the time-averaged energy spectra computed in this 
study and those computed by other authors indicate the difference in averaging 
techniques produces no qualitative difference in the spectra. Only the time-averaged 
spectra are computed by averaging the real and imaginary parts. All other averages are 
computed by taking the magnitude first and then averaging. 

Since the vertical modes are not orthogonal, the energy spectrums as a function 
of vertical modes cannot be precisely compared. However, the meridional energy 
spectrum for each vertical mode is exact since the meridional modes are orthogonal. 
Due to the large number of modes associated with each zonal wavenumber (there are 
nine vertical modes and 37 meridional modes associated with each), it is necessary to 
limit the number of modes to be examined. The choice of modes to be examined 
should be based on the energetics and dynamics of these modes. The energy spectrum 
of these modes should give some indication as to which modes (meridional and 


vertical), if anv, are the most important. 


A. ENERGETICS AND DYNAMICS OF PLANETARY SCALE MODES 

The choice of which vertical modes to examine cannot be based solely on the 
energy spectrum determined in this study because the vertical modes are not completely 
orthogonal, although the larger scale modes are nearly orthogonal. Instead, one must 


also rely on previous studies. Little work has been done in the area of the atmospheric 
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energetics of normal modes, although some related work has been done using empirical 
orthogonal functions. 

The partial basis for the choice of which vertical modes to examine comes from 
Kasahara and Puri (1981) and Ko (1985). Kasahara and Puri (1981) analyzed 
Northern Hemisphere National Meteorological Center (NMC) daily data for the 
‘month of January 1977 using completely orthogonal normal modes. They found that 
the Kinetic energy was the largest for the external mode (2=1) and then it generally 
decreased as a function of vertical mode number, although a substantial amount of 
energy existed in the fourth and fifth vertical modes. They also noted that there was 
more kinetic energy in the smaller scale vertical modes than in other studies that used 
empirical orthogonal functions to represent the vertical structure of the atmosphere. 
An important conclusion of their study is that higher resolution (more vertical modes) 
is needed to clearly resolve the planetary scales. This conclusion was based on an 
analysis of the kinetic energy spectra that showed significant kinetic energy in the 
medium scale vertical modes of planetary-scale waves. 

Ko (1985) analyzed 30 days of perpetual January simulations from the National 
Center for Atmospheric Research (NCAR) Community Climate Model. His results, in 
terms of the kinetic energy spectra, were similar to those of Kasahara and Puri (1981). 
Ko also examined the total and potential energy spectra. He found that a large 
percentage of the total and available potential energy was contained in the medium 
scale vertical modes (£=4-6). Vertical modes 4 and 5 were found to contain the largest 
percentage of total and available potential energy. Those modes also contained about 
the same percentage of kinetic energy as the barotropic and first baroclinic modes. 

The results of Kasahara and Puri (1981) and Ko (1985) indicate that at a 
minimum, the 2 = 1-4 vertical modes need to be examined. Kasahara and Puri (1981) 
presented some Statistics in regards to the energy spectra as a function of meridional 
mode number, but no conclusions were made as to the required meridional resolution 
needed to clearly resolve the planetary scales. 

Because the meridional modes used in this study are orthogonal for a given 
vertical mode, it is possible to examine the energy spectra as a function of meridional 
mode number and draw definite conclusions. The energy spectra for a time average of 
the sum of symmetric (odd indexed) Rossby meridional modes of zonal wavenumber 1, 
2 and 3 for the nine vertical modes of the model are given im Figs. 35.) ="5%53ie 


energy spectra for the large-scale vertical modes (2= 1-3) have a peak in the energy 


around meridional mode number 1-6 with a rapid decrease there-after. The spectra for 
the medium-scale vertical modes (2=4,6) are much flatter than the spectra for the 
large-scale vertical modes. The £2=4 mode still shows a slight peak around meridional 
mode number 9, but the rapid decrease in the spectrum does not occur until meridional 
mode number 19. The spectrum of the 2=5 mode is generally flat until meridional 
mode number 33, while the spectra of the £=6 mode is generally flat through 
meridional mode 57. The general pattern for the £=7 mode is that the energy is 
spread uniformly over most meridional scales. The smallest vertical scale modes 
(2=8,9) have a minimum in energy near meridional mode 17 and then the energy 
increases until mode number 41 where the spectra become relatively flat. The 
corresponding energy spectra for the antisvmmetric (even indexed) Rossby miodes (not 
shown) are very similar. These results are qualitatively consistent with those of Ko 
(1985) and Kasahara and Puri (1982), although there are of course quantitative 
differences. These differences can be attributed to the difference between the modes 
used, the averaging technique used and to the differences in the data analyzed. It 
appears from the energy spectra that no clear choice of which meridional modes to be 
examined can be made. While it is true for certain vertical modes that energy is 
concentrated in particular meridional modes, this pattern does not hold for all vertical 
modes. Clearly some additional criteria are required to reduce the number of modes 
examined. 

It is possible to establish some additional criteria based on the dynamics of the 
different modes. The focus of this studv is on planetarv scale waves and these waves 
have certain dynamical properties that were illustrated by the scale analysis presented 
in Chapter 4. The analysis using the single length scale and the fast time scale (this 
time scale is appropriate for the large-scale barotropic meridional modes, 2=1, n=0-3) 
shows the time tendency to be driven solely by linear effects. The multi-length scale 
analysis using the fast time scale also indicated that the linear terms should be the 
same order of magnitude as the largest term in (3.62). The scale analysis for the slower 
time scales also indicated that the linear term was the same order as the largest term. 
This slow time scale is appropriate for the baroclinic modes. 

The importance of the linear term (same order magnitude as the largest term) can 
be used as a criteria for the selection of the modes to be examined. Only those modes 
that show the linear term to be important should be examined. The amplitudes of the 


terms in (3.62) are given as a function of meridional and vertical niode numbers in 
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Figure 5.1. Energy spectra for an average of zonal wavenumbers 1-3 for (a) the 
barotropic (2=1) rotational mode and (b) Ist baroclinic rotational mode (£=2) as a 
function of the odd (symmetric) meridional mode numbers. 
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Figure 5.2. As in Figure 5.1, except for (a) 2=3 and (b) 2=4. 
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Figure 5.3. As in Figure 5.1, except for (a) 2=5 and (b) 2=6. 
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Figure 5.4. As in Figure 5.1, except for (a) 2=7 and (b) 2=8. 
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Figure 5.5. As in Figure 5.1, except for 2= 9. 
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Figs. 5.6 - 5.10 for a time average of the symmetric zonal wavenumber | case (results 
are similar for wavenumbers 2 and 3). The importance of the linear term 1s evident for 
the n S 3 meridional modes for the first four vertical modes. The linear term is 
important up to meridional mode 8 for the barotropic mode, and up to meridional 
mode 4 for the third vertical mode. For the smaller scale vertical modes (£=7-9) 
modes, the linear term is no longer important for any meridional scale. In general, it 
can be said that the importance of the linear term decreases with increasing vertical 
mode number. Based on the meridional wavenumber arguments given above and the 
analysis of the importance of the linear term as a function of vertical as well as 
meridional mode number it would seem reasonable to limit this study to the 
examination of the first four meridional modes (n= 0-3). 

From Figs. 5.1 - 5.5 it can be seen that the energy contained in then S 3 modes 
decreases as a function of vertical mode. In general the energy contained in then S 3 
modes is of the order 1 - 107! while the energy contained in the n S 3 modes for £ 2 
5 modes is of the order 10°? - 10° The comparison of the amount of energy contained 
in the n S 3 meridional modes is not exact (the vertical modes of the NOGAPS model 
are not completely orthogonal). However, the variation of the energy in these modes 
with vertical mode is very consistent with the results of Kasahara and Puri (1982) and 
Ko (1985). Based on the above comparison of the energy contained in the different 
vertical modes for the n S 3 meridional modes and the consistency of these results 
with those of other authors, it would seem reasonable to limit this study to the 
examination of the first four vertical modes of the n= 0-3 meridional modes. 

Another reason for not examining the smaller scale vertical modes (which cannot 
be addressed by the analysis of Figs. 5.1-5.5 and 5.6-5.10 is the possible dependence of 
the adiabatic nonlinear on the divergent part of the motion. Errico (1984) has shown 
for the long-time solution to a primitive equation model that N(G*R) is the same 
order of magnitude as the N(R*R). The N(G*R) term depends highly on the divergent 
part of the motion. Since the divergent part of the NOGAPS analyses is just the 
model first-guess divergence, any results obtained concerning nonlinear interactions 
may be more indicative of the model rather than the actual atmosphere. The 
dependency of the interactions on the model is an important subject in its own right, 
but it is a subject beyond the scope of this work. 

The choice of of the modes to be examined (the 2=1-4 , n=0-3 modes) does 


eliminate some modes with energies similar to the chosen modes and modes that show 
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the linear term to be important, but the set chosen should be sufficient for the purpose 
of this study which is to show that nonlinear interactions with cyclone waves can be an 
important factor in the dynamics of planetary waves. 

Before examining the results of the experiment that is designed to determine the 
importance of synoptic-scale nonlinear interactions, it is necessarv to establish the 
importance of the adiabatic nonlinear term in the dynamics of planetary waves. While 
it may be shown that synoptic-scale interactions are an important part of the adiabatic 
nonlinear term, this effect would not be significant unless this term 1s important in 
general (of the same magnitude as the other terms). As can be seen from Figs. 5.6 - 
5.10 the adiabatic nonlinear term is at least the same order as the largest term for a 
given meridional mode and it is often the largest term. For the largest scale vertical 
modes, all of the terms (linear, nonlinear adiabatic, diabatic and time tendency) are 
nearly the same magnitude for the largest scale meridional modes (n= 1-3). This result 
iS interesting considering the single length scale analysis carried out using the fast time 
scale (which is appropriate for the 2=1, n=0-3 modes) showed the time tendency to 
be driven by the linear term. The multi-scale analysis did show that the nonlinear term 
could be important. For most of the remaining meridional modes, the time tendency 
for the largest vertical scales is driven by the adiabatic nonlinear term, with the diabatic 
term becoming important for the smallest scale meridional modes. 

The importance of the diabatic term for the largest vertical scale meridional 
modes has not been shown before. In the only other study of this kind, Errico (1984) 
did not investigate the planetary scales in detail and his focus was mainly on the 
balance for the gravity modes. In addition, he averaged modes by frequency, which 
could have obscured some of the detail of the planetary scales. The results that he 
presented did show the adiabatic nonlinear term to be important for the higher 
frequency Rossby modes, but he did not show the diabatic term to be important for 
the Rossby modes except for the higher vertical modes. Although it beyond the scope 
of this study to examine in detail the affects of diabatic terms on planetary waves, the 
importance of the diabatic term, as shown in Figs. 5.6 - 5.10, does deserve some 
comment. 

It is not surprising that this term might be important for planetary scales given 
the general belief that the quasi-stationarv component of planetary waves is forced in 
part by differential heating due to land-sea contrasts. A note of caution must be added 


here. As is the case with most models, the diabatic processes of the NOGAPS model 
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Figure 5.6. Magnitude of AC /At (dotted line), -imC, (solid), N (dash-dot), and Q. 
(dashed) for (a) the barotropic (€=1) rotational mode and (b) the Ist baroclinic 
Rotational mode (£=2) as a function of meridional mode number (n) for zonal 
wavenumber lI. 
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Figure 5.7. As in Figure 5.6, except for (a) 2=3 and (b) 2=4. 
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Figure 5.8. As in Figure 5.6, except for (a) 2=5 and (b) £=6. 
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As in Figure 5.6, except for (a) 2=7 and (b) 2=8. 
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Figure 5.10. As in Figure 5.6, except for 2=9. 
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are parameterizations of the actual atmospheric diabatic processes. Thus results 
obtained here may be different from other models with different parameterizations. In 
addition, the diabatic term is not decomposed into different components (sensible 
heating, latent heating, frictional dissipation, radiational heating), so it is not possible 
to determine which diabatic processes contribute the most to this term. For the larger 
meridional scales, it is possible that the majority of this term is due to heating. For the 


smaller scales, dissipation may be important. 


B. THE EFFECT OF SYNOPTIC SCALES ON PLANETARY SCALES 

Based on the dynamics and energetics of the different modes and on scale 
considerations, the focus will be on zonal wavenumbers 1-3, vertical modes 1-4 and 
meridional modes 0-3. For this part of the experiment, the data (u. v, T, In p.) for 
each of the 19 davs of LLU TC NOGAPS analvses taken every 5 davs from 16 Janta 
1986 to 16 April 1986 are spectrallv filtered. The data are taken onlv everv 5 days so 
that each day 1s statistically independent. Also, graphs of the data taken everv 5 days 
are easier to interpret than those that have daily data plotted. When average 
magnitudes and energies are computed later, 85 days of data taken over the same 
period will be used. 

To eliminate the synoptic waves. the data are filtered by transforming to spectral 
space using a fast Fourier transform routine, setting the coefficients of zonal 
wavenumbers 7-15 to zero and transforming back to physical space. The magnitudes 
of the adiabatic nonlinear term (N,,) from both the unfiltered and filtered data are 
computed and the difference between the nonlinear terms from the unfiltered and 
filtered data is calculated. The magnitude of this difference represents the contribution 
of the synoptic scales to the nonlinear term of the given planetary scale, hereafter 
referred to as Nj. 


Figures 5.11 - 5.14 are plots of N, and N,, versus time for an average of 


S 
meridional modes n=(-3 of the £= 1-4 modes of zonal wavenumbers m=1-3. The 
general pattern that emerges from an examination of Figs. 5.1] - 5.14 is that for all 
vertical modes N,, 1s the same order of magnitude as N_. but it is generally smalter 
than N,, (here after referred to as the total nonlinear term). On no day is N,,, larger 
than N,, (which could occur if the interactions from waves other than m= 7-15 acted to 
oppose the interactions from the m= 7-15 waves). On a number of days, the synoptic- 
scale contribution to the total nonlinear term is significant (N 


‘ns 1S a large percentage 
cian): 
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For a given day, the significance of the synoptic-scale contribution to the total 
term can be illustrated by the ratio of N,, to N,. Figures 5.15 - 5.18 are plots of the 
percent magnitude of the synoptic-scale contribution to the total nonlinear term as a 
function of time. The synoptic-scale contribution to the magnitude of the total 
nonlinear term will be considered significant for any given dav if this contribution is at 
least 30 percent of the total. This is somewhat arbitrary, because a smaller value could 
be important if AC/At is smaller than the nonlinear term. However, the 30% value 
does highlight days on which the synoptic-scale contribution to the nonlinear term 1s 
large. The number of days where the magnitude of the synoptic-scale contributions 1s 
greater than 30 percent varies with vertical mode. The number of significant days 1s 
greatest for fourth vertical mode. In general, the larger scale vertical modes (2= 1-3) 
have the smallest number of significant davs. The strength of the svnoptic-scale 
contribution to the total nonlinear term varies with vertical mode number in the same 
manner as the number of significant days. The synoptic-scale contribution to the total 
is strongest for medium-scale fourth vertical mode. For this mode, the svnoptic-scale 
contribution to the total nonlinear term exceeds 50 percent for a number of days. 
There is a slight variation of the synoptic-scale contribution to the total nonlinear term 
with zonal wavenumber. Zonal wavenumber 3 1s the most affected zonal wave 
number. The fourth vertical mode for this zonal wavenumber exceeds the 30 percent 
criteria for 12 of the 19 days and on three of these days the svnoptic-scale 
contributions exceeds 50 percent. 

From the above analvsis, it is clear that synoptic-scale interactions can be 
significant for a particular day, but the analysis does not show the time average effect 
of the synoptic scales on planetary scales. Table 6 gives the average magnitude of the 
difference between the total nonlinear term and the nonlinear term computed from the 
filtered data as a percentage of the average magnitude of the total nonlinear term. The 


averages for these terms are computed as 
es 7 
Nae =!— Ny lth )N; (t 
re d / Nn (te) Ne (te) 


=lilna= 


(5.1) 

Where K is the number of davs in the data set and four is the number of meridional 
modes. To insure greater statistical significance, a larger sample size is used to 
compute the averages given in Table 6. For these averages, the total term (N,) and 
the synoptic part (N,,) are determined from 85 days of 12Z initialized NOGAPS 
analyses taken every day from 16 January 1986 to 16 April 1986. The following days 


were missing from the data set: 17, 21 and 23 January and 9 and 21 March. 
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Figure 5.11. Magnitudes of N_ (solid) and the synoptic scale contribution to N 
(dashed) as a function of time for zonal wavenumbers (a) I, (b) 2, and (c) 3 for the 
barotrophic mode (£= 1). 
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Figure 5.12. As in Figure 5.11, except for 2=2. 
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Figure 5.13. As in Figure 5.11, except for 2=3. 
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Figure 5.14. As in Figure 5.11, except for 2£= 4. 
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Figure 5.15. Synoptic-scale contribution to the total nonlinear term for the barotropic 
mode (£=1) as a percentage of the total nonlinear term for zonal wavenumbers 
m= [-3 (a-c). The 30% level is indicated for use as a benchmark. 
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Figure 5.16. As in Figure 5.15, except for 2= 2. 
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Figure 5.17. As in Figure 5.15, except for 2= 3. 
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Figure 5.18. As in Figure 5.15, except for £=4. 
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Table 6. Time averaged (over 85 days, taken between 16 January and 16 April 
1986) synoptic-scale contribution to the total nonlinear term as a percentage of the 
total for an average of meridional modes n=0-3. 


Vertical Zonal Wavenumber 
mode index 
a4 2 3 
i 20 18 ey 
Z Ly 14 16 
3 24 25 24 
= 35 28 42 


Fora given zonal wavenumber, the magnitude of the time averaged svnoptic- 
scale contributions is a largest percentage of the total nonlinear term for the fourth 
vertical mode. The nonlinear term of the second vertical mode appears to have the 
smallest contribution from synoptic scales. For the medium-scale vertical modes, zonal 
wavenumber 3 has the largest svnoptic-scale contribution to the planetary scales in 
terms of the percent magnitudes. 

It is clear from the above analvsis that the svnoptic-scale contributions to the 
magnitude of the total nonlinear term can be significant for a given day (N_, can be as 
large as 60% of N,,) and are significant in a ume-averaged sense. 

The result that the svnoptic scales tend to have a significant mmpact on the 
magnitude of the total nonlinear term 1s important, but it 1s not the whole story. For 
the most affected mode (2=4, m= 3), the nonlinear interaction from synoptic scales is 
42%. of the total. However, this leaves nearly 60° of the total unaccounted for. The 
rest of the term 1s due to interactions with scales other than those represented bv 
wavenumbers 7-15. This includes interaction of the planetarv scaies with the zonal 
mean state (m=Q). Although these other interactions are significant, it 1s bevond the 
scope of this studv to examine them. The purpose of this studv 1s to establish that 
interactions with cvclone waves can be an important factor in the dynamics of 
planetary waves. 

Another point that has not been considered is the relationship of the davs when 


Nps 1S a large percentage of N,, to the general importance of N,, on those days. It is 
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important to establish that on the days when the synoptic-scale contribution to the 
total nonlinear term is a large, the total nonlinear term is of importance, Le., it is of the 
Same magnitude as the other terms in (3.62). To accomplish the above, one must 
examine the non-time averaged variations of the terms given in (3.62). [n addition, 
such an examination will also serve to illustrate the dvnamical differences between the 
different vertical modes at these planetary scales. 

The terms in (3.62) are plotted in Figs. 5.19 - 5.22 for an average of meridional 
modes n=0-3 versus time for the data taken every five days. The time variations of 
barotropic mode (Fig. 5.19) are quite interesting. For all three zonal wavenumbers, the 
general dynamic pattern is a tendency for balance between the adiabatic nonlinear term 
and the linear term. The time tendency of these modes is generally less than the 
nonlinear or linear term. and the diabatic term also has a small contribution to the 
barotropic mode. On the davs on which the nonlinear and linear terms are not in 
partial balance, the time tendency is driven by the linear term. For these days, it 
appears as if the heating is partially balancing the nonlinear term. For all the days 
examined, the contribution of the nonlinear to the barotropic mode term 1s at least the 
Same order of magnitude as the other terms. It was previously shown (Fig. 5.15) that 
the synoptic-scale contribution to the total nonlinear term for the barotropic mode is 
small. However, the synoptic-scale contribution to the overall dynamics of the 
planetary-scale barotropic modes could be more important than the small magnitude 
might indicate. This contribution could be an important factor in maintaining the 
partial balance found in this mode or it could be the dynamic ingredient that keeps 
these modes from obtaining a complete balance. The above could also be true of the 
diabatic term which also has a small magnitude for the barotropic modes. 

The above description of the dynamical nature of the barotropic mode is also 
generally true for the second vertical mode (Fig. 5.20). The dynamics of the third (Fig. 
5.21) and fourth (5.22) vertical modes are quite complex. No one term or two terms 
seem to dominate. The diabatic term is of increased importance for these modes, but it 
is by no means dominant. The linear term is of lesser importance for these modes than 
it is for the first two vertical modes, but it still makes a significant contribution to the 
time tendency. There is no hint of a simple balance between two leading terms for 
these modes. For both of these modes, it generally true that the diabatic, adiabatic 
nonlinear and linear terms are additive so that the time tendency is the largest term. 


Two important points are illustrated by the above analysis 
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Figure 5.19. Magnitude of AC /At (dotted), -i@C_ (solid), N (dash-dot), and Q. 
(dashed) for the barotropic (£= 1) modes of zonal wavenumbers 1-3 (a-c) for each fifth 
day of the period. 
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Figure 5.20. As in Figure 5.19, except for £= 2. 
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Figure 5.21. As in Figure 5.19, except for 2= 3. 
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Figure 5.22. As in Figure 5.19, except for 2=4. 
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e For the vertical mode that 1s most affected by svnoptic-scale interactions 
(£=4), the dynamics are the most complex. The nonlinear term is on average 
important for this mode, but this importance varies significantly over the period 
being examined. 


e The time tendencies of the barotropic and first baroclinic modes are generally 
smaller than the nonlinear term due to the partial balance between the linear 
and nonlinear terms. Thus, synoptic-scale contributions can be more important 
in determining the time tendency of these modes than their small percentage of 
the total nonlinear term may indicate. 


To establish the importance of synoptic-scale interactions for these modes, it ts 
necessary to compare the magnitude of the synoptic-scale contribution to the total 
nonlinear term to the magnitude of the time tendency. The magnitude of the synoptic- 
scale contribution to the total nonlinear term as a percentage of the magnitude of the 
time tendency are given in Figs. 5.23 - 5.26 for the £=1-4 vertical modes of zonal 
wavenumbers 1-3. In addition, the synoptic-scale contribution to the total nonlinear 
term is plotted to allow comparisons on the days when the synoptic scales make up a 
large percentage of the total nonlinear term. For the barotropic mode (and especially 
for the barotropic mode of zonal wavenumber 3), the synoptic-scale contribution is a 
larger percentage of the time tendency than it is of the total nonlinear term. For 
example, the svnoptic-scale contribution to the total nonlinear term of the barotropic 
mode of zonal wavenumber 2 on the 26th and 31st of January is only about 12-14% of 
the total nonlinear term, but 38-40% of the time tendency. This result implies that for 
the barotropic modes, the synoptic-scale contributions are more important than the 
small value of N..,N,, would indicate. The result is also true, although to a lesser 
extent, for the second vertical mode (Fig. 5.24) of zonal wavenumbers | and 2. For 
vertical modes 3 and 4 (Figs. 5.25 - 5.26), the general pattern 1s that the synoptic-scale 
contributions to the total nonlinear term are a larger percentage of that term than of 
the time tendency. This is an indication that the nonlinear term 1s no longer the 
dominant term for these modes. Rather, the nonlinear term it 1s only one of three 
terms that contribute to the time tendencv and that these terms are additive. The lack 
of any balance for the £2=3,4 modes is an indication that the linear term is no longer 
large enough to balance the nonlinear term. However, the fourth vertical mode is still 
the most affected by synoptic scales. The number of days for which the magnitude of 
the synoptic-scale contributions to the total nonlinear term is greater than 30 percent 
of the magnitude of the time tendency is generally less than the number of days for 


which the synoptic-scales contributions are greater than 30 percent of the total 
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nonlinear term. However, there are still more significant days for the fourth vertical 


mode than anv other vertical mode. 


C. ENERGY RELATIONSHIPS 

It is clear from the above analysis that synoptic scales can play an important role 
in the dynamics of planetary scale waves. It is also clear that these interactions are 
important in a time mean sense (Table 6). The analysis in the previous section is very 
good at highlighting this importance, but it does not show how the energy flows 
between synoptic scales and planetary scales. Are the scale interactions acting to 
decrease or increase the amplitudes (energy) of these modes? This question will be 
addressed in the this section. It will also be very important to determine the average 
energy transfers that can be compared with other studies (e.g. Tanaka et al., 1980). 
Most of these studies considered either kinetic or potential energy transfers, while total 
energy energy transfers are treated in this study. 

The sum of kinetic plus available potential energy in a particular mode is given 
simply by 

Cs (S22) 

where C,, is the amplitude of the mode designated by the index n and ( )" indicates the 
complex conjugate. Thus, an equation for the time tendency of the total energy of a 
given mode can be obtained by multiplying (3.62) by C,,* , and adding this to the 


product of the complex conjugate of (3.62) times C,. The contribution of N,, to the 


time tendency of the energy of a mode n for m # 0 is given by 
2Re(C,N_) (5.3) 


where Re indicates the real part of the expression. The relation given in (5.3) indicates 
that the energy tendency generated by N., or N,,, depends on the phasing between e 
and N, or N,.. If the synoptic-scale forcing is out of phase with the planetary scale 
mode, then a large magnitude of N,, does not necessarily imply a large energy 
transfer to the planetary mode. However, this might still indicate an important 
contribution to the phase speed of the planetary wave. 

The synoptic-scale contribution to the energy tendency produced bv the 


nonlinear term is computed by taking the difference of the tendencies as computed in 
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Figure 5.23. Synoptic-scale contribution to the total nonlinear term of the barotropic | 
mode (£= 1) as a percentage of the magnitude of the time tendency (solid) and of the 
total nonlinear term (dashed) for zonal wavenumbers m= I-3 (a-c). 
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Figure 5.24. As in Figure 5.23, except for 2= 2. 
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Figure 5.25. As in Figure 5.23, except for £= 3. 
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Figure 5.26. As in Figure 5.23, except for 2=4. 
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(5.3) from the filtered and unfiltered data. The energy tendency that is associated with 
the synoptic-scale ccntribution to the nonlinear term is shown in Figs. 5.27 - 5.30 for 
the 2= 1-4 vertical modes of zonal wavenumbers 1-3. Also included on these plots are 
the contributions of the total nonlinear term. For the barotropic mode of zonal 
wavenumber I (Fig. 5.27), there is clearly a bias for the synoptic-scale contribution to 
the energy tendency to be positive. That is, the synoptic scales are interacting such 
that they are increasing the energy of the barotropic mode on 16 of the 19 days that 
were analyzed. The above pattern is evident to a lesser extent for the barotropic 
modes of zonal wavenumbers 2 and 3. For the 2=2,3 modes (Figs. 5.28-5.29), there is 
a tendency for the synoptic scales to take energy awav from the planetary scale. This 
tendency appears to be strongest for the 2=3 modes. For example, on only two of the 
19 davs are the svnoptic scales transfering energy to zonal wavenumbers 3 and 2 scales. 
The fourth (Fig. 5.50) vertical mode has the largest contribution from synoptic scales 
to planetary scales in terms of the percentage of the total energy tendency due to the 
nonlinear term. For most of the days examined, the synoptic-scale contribution to the 
energy tendency is positive for all three zonal wavenumbers. 

Table 7 gives the time-averaged synoptic-scale contribution to the energy 
tendencies of the nine vertical modes of zonal wavenumbers 1-3. Energy flow from 
synoptic scales to planetary scales would be indicated by a positive value. Once again, 
tO insure greater statistical significance the averages in Table 7 are computed from 85 
days of 12UTC initialized NOGAPS analyses 

The general pattern that is evident from Table 7 is that barotropic (2=1), and 
third baroclinic (2=4) modes have a positive contribution from synoptic scales to the 
energy tendency for all three zonal wavenumbers. The synoptic scale contribution for 
the £=2,3 modes are consistently negative for all three zonal wavenumbers. The 
svnoptic-scale contribution to the energy tendency is generally largest for the 
barotropic mode and smallest for the 2=4 mode. The relative importance of the 
energy flow from synoptic to planetary scales for a particular vertical mode is not 
truely indicated by energy tendencies given in Table (7). For example, the barotropic 
mode has a larger synoptic-scale contribution to energy tendency than does the 2=4 
mode, but the barotropic mode contains a larger amount of energy than does the 2=4 
mode. However, the most important point to be made from Table (7) is not the 
strength of the synoptic-scale contribution, but rather the pattern of this contribution. 
Based on the the pattern of energy flow from synoptic to planetary scales, one can 


hypothesize some possible mechanisms for these interactions. 
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Figure 5.27. Energy tendency due to the total nonlinear term (solid) for the barotropic 
mode (£=1) and the synoptic scale contribution to the energy tendency generated by 
the nonlinear term (dashed) for zonal wavenumbers m= 1-3 (a-c). 
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Figure 5.28. As in Figure 5.27, except for 2= 2. 
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Figure 5.29. As in Figure 5.27, except for 2=3. 
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Figure 5.30. As in Figure 5.27, except for 2=4. 
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Table 7. Time averaged (over 85 days, taken between 16 January and 16 April 
1986) svnoptic-scale contribution to the energy tendency (J Kg! s"!) generated by 
the total nonlinear term for an average of meridional modes n=0-3. 
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The analvsis of the energy flows in this study between synoptic scales and 
planetary scales is unique in two wavs: The analvsis 1s done considering the total 
energy flow and it 1s specific to interactions between svnoptic waves and planetary 
waves. Tanaka et al. (1986) examined kinetic energy flow for groups of modes, but 
they did not look at the total (kinetic and potential) energy flow and their emphasis 
was on zonal mean-eddy interaction. A important result of their study was that the 
zonal mean barotropic mode gained energy from the baroclinic modes of higher 
wavenumbers. Tanaka et al. indicated that the kinetic energy source for the higher 
wavenumber baroclinic modes came from the conversion of potential energy via 
baroclinic instability. 

The results given in Table 7 show that the barotropic modes of all three zonal 
wavenumbers gain total energy from synoptic scales while the planetarv-scale 
baroclinic modes (£= 2,3) were losing total energy to synopttc scales. These results are 
simular to the results of Tanaka et al. (1986) if one thinks of a local basic state having a 
projection on planetary scales. Here I use the term basic state to refer to the state 
obtained by zonallv averaging over a restricted domain, such as the wavelength of one 
cyclone. The above is in contrast to a zonal basic state that implies averaging around 
an entire latitude circle. Because of the large scales of planetary waves, a local basic 
state may have a large projection on the planetary scales. Thus, the synoptic scale 
£=2 and 3 vertical modes may be gaining potential energy from this basic state, 


converting this potential energy to kinetic energy via baroclinic instability and then 
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transfering this kinetic energy back to the barotropic (2=0) component of basic state. 
In addition, the synoptic scale 2=2 and 3 vertical modes could at the same be 
transfering energy to the 2=4 planetary scale modes. 

The hypothesized mechanism for the energy flow described above is also 
consistent with the results obtained from the analytic model presented in Chapter 4. 
For example, if the zonal wavenumber 6 mode for the first baroclinic mode (2=1 in 
the notation of the simple analytic model, £2=2 for the NOGAPS analysis) were to 
interact with the zonal wavenumber 7 mode for the second baroclinic mode (2=2 for 
the analytic model, and 2=3 for the NOGAPS model) then the zonal wavenumber 1 
mode for the first baroclinic mode (2=1 for the analytic model, 2 for the NOGAPS 
model) and third baroclinic (2=3 for the analvtic model, 2=4 for the NOGAPS 
model) modes would be affected. Also, the first two baroclinic modes could interact 
with themselves to affect the barotropic mode. 

Since the results presented in Table 7 are for the total energy flow it is not 
possible tell how much of the energy flow is kinetic energy and how much is potential 
energy. The results from the simple analytic model presented in Chapter IV, indicate 
that the majority (if not all) of the energy flow to the barotropic mode should be due 
to Kinetic energy transfer while both kinetic and potential energy transfer are possible 
for the baroclinic modes. In the following section the relative importance of energy 
transfer through the momentum advection terms (kinetic energy transfer) to the energy 


transfer through the mass advection terms (potential energy transfer) will be examined. 


D. SOME MECHANISMS OF SYNOPTIC-SCALE INTERACTIONS 

It is clear from the analysis of the previous sections that synoptic scales can have 
a significant impact on the dynamics and energetics of planetary scales. What is not 
clear are the mechanisms by which the synoptic scales are interacting with the 
planetary scales. Are mass field interactions (e.g. temperature advection) the primary 
mechanism (Gall et al., 1979) or are momentum field interactions (e.g. momentum) 
advection more important? One of the reasons why the mechanisms of the interaction 
are not clear is that the adiabatic nonlinear term is a combination of of terms from the 
momentum equations and the thermodynamic energy equation. That is, the nonlinear 
term contains the effects of temperature advection as well as momentum advection. 
What is needed is a way to separate these effects. The method used to separate these 
effects is simple and direct. That part of the synoptic-scale contribution due to 
momentum field interactions is isolated in the following way: 


°e. The data are filtered as described in the beginning of this Chapter. 


94 


e The nonlinear term for both the filtered and unfiltered data is computed by first 
computing the discrete time tendency of the coefficients via an adiabatic one- 
time step integration of the NOGAPS model, and then subtracting the linear 
term (1@C) from it. However, in this step the term representing the horizontal 
momentum advections is deleted. 


e The energy tendency due to the nonlinear term with the momentum advections 
deleted is computed using (5.3) from both the unfiltered and filtered data. 


e The energy tendency due to synoptic scales with and without the momentum 
advections is determined by subtracting the energy tendencv computed from the 
filtered data set from the energy tendency computed from the unfiltered data for 
each case. 


e The energy tendency due to synoptic scales that comes through the momentum 
advections is computed by subtracting the energy tendency due to svnoptic 
scales computed without momentum advections from the svynoptic-scale 
contribution to the energy tendency computed from integrations made with all 
of the advection terms included. 


Care must be taken in using the above method. Since the svnopuc scale 
contribution to the nonlinear term can often be a small difference between two large 
terms (the nonlinear terms computed from the filtered and unfiltered data) care must 
be taken to insure that the deletion of the momentum advection terms does not 
produce a large increase in time tendencies. A large increase in the time tendencies 
may mean that one would be trving to determine a value as a very small difference 
between two very large numbers. If this difference is smaller than the accuracy of the 
two values being subtracted, then the difference 1s not reliable. The error produced bv 
the filtering process is such that only 3 to 4 significant digits are maintained after 
filtering. Thus if the original difference is small (an order of magnitude smaller than 
the original terms), and if the deletion of the momentum terms increases the time 
tendency by two orders of magnitude or more then the results would be questionable. 

The other potential source of error in this method arises because the model 
equations are in a flux form so that deleting the term that represents the momentum 
advections also means deleting other terms (which enter via the continuity equation) 
which do not involve momentum advections. The continuous flux form of the 
momentum advection term in the sigma coordinate of the NOGAPS model is 

envy ) + a“ (5.4) 
Where 7 is defined by (3.10). The discrete form of the first term in (5.4) is deleted from 


the NOGAPS model one-time step integrations as part of the analysis procedure which 
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is designed to isolate interactions through the momentum advection term. The first 


term in (5.4) can be written 


mV -VV)+2VVV +V(V-Vzx). 


The variable % may be written as 


= t 
Tat (5.6) 
where 7 is a constant over the globe and # is a small deviation from this constant. In 


general, 7>>7' except near high mountains. Thus the first term in (5.4) 1s 


approximated by 


mV - VV). 


Cay 


( 


The third term in (5.5) will in general be much smaller than (5.7) except near 


.7) 
mountains. The second term in (5.4) and (5.5) will be small for the modes being 
considered because they represent interactions between the divergent and rotational 
parts of the motion (Errico, 1984). 

Based on the above arguments concerning the smallness of the neglected terms, it 
can be seen that most of the nonlinear interactions that may occur through the flux 
form of the momentum advection term will be due to (5.7). 

This method for determining the relative importantance of momentum and mass 
field interaction 1S most accurate when the synoptic-scale contribution to the total 
nonlinear term is large and when the deletions of the advection terms do not produce 
extremelv large changes (at least two orders of magnitude) in the time tendencies of the 
modes. Deletion of the momentum advection terms does not produce extremely large 
changes to the time tendency of the modes. The above is not true when the 
temperature advection terms were deleted. Deletion of these terms produces large 
changes in the time tendencies of most of the baroclinic modes. Another point to 
consider when deleting nonlinear terms is that one must not delete terms that involve 
the mean state that the primitive equations are linearized about. For example. deietion 


of the vertical advection of temperature would delete the linear basic term 


Go" (5.8) 
Such a deletion would mean that the frequencies that are determined as eigenvalues of 


the linearized equations, would no longer be valid. As result, the linear term 10C 
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would not be accurate. Since the linearization about a mean state of rest contains no 
horizontal advections, the deletions of the horizontal momentum advections will not 
alter the frequencies. Because deletion of the temperature advection produces large 
changes in the time tendency of the modes, the interactions that are due to mass field 
interactions are assumed to be the values computed with the momentum advection 
‘turned off. These values do not represent just the horizontal temperature interactions 
since in sigma coordinates there are a number of nonlinear mass interaction terms e.g. 
the advection of surface pressure and the horizonal gradient of surface pressure. 
However, the temperature advection would be expected to comprise a substantial part 
of these nonlinear terms. 

The synoptic-scale contribution to the energy tendency generated by the 
nonlinear term computed with all advections, and the difference between the synoptic 
scale contribution with all advection and that quantity computed without momentum 
advections are plotted in Figs. 5.31 - 5.34. The difference represents the portion of the 
synoptic scale contribution to the energy tendency that is due to the momentum 
advections. For the barotropic mode (Fig. 5.31), the synoptic scale contribution to the 
energy tendency generated by the nonlinear term is mainly due to momentum 
interactions. This is especially true when the synoptic scale contribution to the energy 
tendency is large. For example, compare the values of the total synoptic scale 
contribution to the energy tendency with that part due to the momentum advection 
terms on: 26 January, 17 and 27 March for zonal wavenumber one; 17 March and 6 
Apnil for zonal wavenumber two; and 26 January, 27 March and | Apmril for zonal 
wavenumber three. Momentum advections are a slightly more important mechanism 
than mass field interactions for the second (Fig. 5.32) and third (Fig 5.33) vertical 
modes. However, the mass field interactions can be the larger contribution for some 
days ( 20 February for 2=3, m=1; 26 January, 25 February, and 27 March for 2=2, 
m=2; and 5 February for 2=2, m=3). Mass field interactions are generally the 
dominant mechanism for the fourth (Fig. 5.34). This is especially true for days with a 
large energy tendency due to synoptic-scale interactions. However, there are also a a 
few days where momentum advections make a substantial contribution. 

The results given in Table 7 show that the barotropic mode of all three zonal 
wavenumbers gains total energy from the synoptic scales. The results presented in Fig. 
(5.31) indicate that this total energy transfer is mainly due to momentum interactions. 


Thus, the energy being transfered from synoptic to planetary scales is kinetic energy. 
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Figure 5.31. Synoptic scale contribution to the energy tendency generated by the 
nonlinear term for the barotropic mode (£=1) with all advections (solid) and the 
contribution due to synoptic scale momentum advections (dashed) for zonal 
wavenumbers m= 1-3 (a-c). 
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Figure 5.32. As in Figure 5.31, except for 2=2. 
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Figure 5.33. As in Figure 5.31, except for £=3. 


100 


-7 
*%10 
15.0 
ae | 











Sp 
we 
iS: 
= 
> 
¢3 © 
>. 
Oo) Ww 
=a) 
roll 
Ge} 
nc 
Renn = 
Ey) 
A y 
Fa 
es 
16 21 2631 S$ 19152025 2 7121722277 1 6 11 16 
JAN FEB MAR APR 
1986 
DAY 
> ° 
Ono Oo 
x 
x ' eS 
Ny 
3 Dp 
Sa wo 
™~ ac ~ = 
- ? ‘ 
ry 
v4 
> > ot 
ey. © (=) De ; 
Par = N iN aa ‘ / 4 
| Ge ‘\ POEs ’ 
=) CS i \ f \ 
<— ae , ; ‘i 
Oj CJ ; “hoot! ee ‘ 
E39 Fo ls i Ye 
~~ 2 Ol 
C5 oe ‘oy Ca \/ 
9 al re a) 
C13 ve al 
a =z° nN 
Lin ooo Hal 
] i 
Herein S| S 1915 202 2 7 12 17 22 27 : 6 11 16 “Stee 5 10152025 2 7 12172227 1 6 1116 
JAN FEB MAR APR SAN PES MAR APR 
1986 1986 
DAY DAY 


Figure 5.34. As in Figure 5.31, except for 2=4. 
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This result is consistent with the hypothesis presented in the previous section which: is 
that the relation of the svnoptic-scale waves to the planetaryv-scale waves can be 
likened to the relationship of an eddy to local mean flow. Assunung this hypothesis 1s 
correct, the transfer of kinetic energy from the svnoptic-scale modes to the barotropic 
planetary-scale modes mav be associated with the process of barotropic stability. 

The results presented above for the barotropic mode are also consistent with 
those of the simple analytic model of Chapter IV. The results from this model indicate 
that mass field interactions would not be possible for the barotropic mode. This is not 
strictly the true for NOGAPS model, but it does appear from Fig. 5.31 that the 
momentum advections are also dominant in the NOGAPS barotropic mode. 

The results presented in Table 7 indicate that on the average the synoptic scales 
are taking energy from the planetary-scale 2= 2 and 3 modes while giving energy up to 
the planetary-scale 2=4 mode. From Figs 5.32-5.34 it is clear that the total energy 
transfer between these scales is due to both kinetic and potential energy transfer. The 
results from the simple analytic model indicated that both types of energy transfer are 
possible for baroclinic modes. Also the results for the £=2-4 modes are not 
inconsistent with the idea that the relationship of the synoptic-scale modes to 
planetary-scale modes can be likened to the relationship of an eddy to a local mean 
flow. There are three possibilities for kinetic and potential energy flow between 
synoptic- and planetary-scale baroclinic modes: 


e Synoptic-scale baroclinic modes can gain both kinetic and potential energy from 
planetary-scale baroclinic modes. 


¢ Synoptic-scale waves can gain potential energy from planetary-scale baroclinic 
modes while losing kinetic energy to those same modes. 


e Synoptic-scale modes can lose potential energy to planetary-scale modes while 
gaining kinetic energy from those same modes. 


The final consequence of the above analogy, is that the type, sign and amount of 
energy transfer from synoptic-scale baroclinic waves to planetary-scale baroclinic waves 
woulu depend in a crucial way on the tilt of the synoptic scale waves’ phase with 
respect to the vertical and horizontal wind shear of the planetary waves. 

The potential energy flow from synoptic scales to the planetary scales of the £=4 
mode is similar to the flow proposed by Gall (1979). Gall proposed that the planetary 
scale waves were forced mainly bv planetary scale variations in the meridional heat flux 
convergence of higher wavenumber modes 1.e., the interaction between cyclone scales 


waves and the local basic state increased the amplitude of the planetary scale 
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temperature perturbation. This potential energy was then converted to kinetic energy 
due to a positive correlation between planetary scale upward motion and temperature. 
It is clear from the results presented in Fig. 5.34 that potential energy is being 
transfered to the £=4 mode of the planetary scale waves. The above result is also 
consistent with the results obtained from the analytic model presented in Chapter IV. 
However, the results of this study also appear to conflict with the results of Gall et al. 
(1979) in that not all of the planetary scale modes (the 2=1 mode for example) are 
being forced mass field interactions. The apparent conflict between this study and that 
of Gall can be explained by the differences in the data rather than the mechanism. In 
Gall’s experiment, synoptic-scale perturbations were allowed to grow from a Zonal 
mean state. Gall argued that in developing waves/perturbations that 
Wave/perturbation velocity (V’) and the wave/perturbation vorticitv (¢’) were 
uncorrelated while V’ and the wave/perturbation temperature were correlated. [n this 
study, the atmosphere data contained a variety of fully developed cyclones. For fully 
developed cyclones it cannot be said that V’G’ are uncorrelated. In fact, general 
circulation theory would suggest that V’C’ would be correlated. In a more recent study 
that was very similar to that done by Gall, Young and Villere (1985) showed that direct 
transfer of kinetic energy from intermediate scales to planetary scales was of equal 
importance to the transfer of potential energy. A possible reasan for this conflict 
between these two similar studies is that the zonal mean state specified by Young and 
Villere was such that they obtained higher growth rates than Gall and their 
disturbances developed faster. I believe that the disturbances in the Young and Villere 
studv developed to the stage where the correlation V’C’ became significant. While it 
may be true for the simple state specified by Gall that the synopticscale forcing of 
planetary scale waves is mainly through the temperature advection term, this does not 
appear to be true for an atmosphere that contains fully developed cyclones. This is not 
to say that the mechanism proposed by Gall is invalid, as this mechanism may be 
active in the atmosphere given the proper distribution and variation of intensity of 
cyclones. However, the direct transfer of kinetic energy could be taking place that may 
Or may not be in the same sense as the potential energy transfer. 

In summary, the analysis of the results presented in this chapter have established 
the following: 


e Synoptic-scale interactions can have a significant impact on the dynamics of 
planetary scale modes. 
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¢ For the barotropic mode (2=1), the main balance is between the linear (Rossby 
wave) terms and adiabatic advection terms, so that the time derivative is smaller 
than either term. In this situation. the svnoptic-scale interactions can have an 
important influence on the planetary scale waves even when they are not a large 
percentage of the total advection term. 


e Synoptic-scale interactions tend to decrease the energy of the planetary 2=2,3 
vertical modes while increasing the energy of the barotropic and 2=4 modes. 


e The synoptic-scale interactions are mainly through the momentum terms for the 
barotropic mode and through both the momentum and mass advection terms 
for the baroclinic modes (£=2-4 modes). The mass field interactions are 
generally dominate for the 2=4 modes. 


There are some important implications of the above results to the forecastability 
of planetary-scale waves. One can easily see that if the strengths and/or phases of 
svnoptic-scale waves are not not forecast correctly then the interactions between 
synoptic and planetary scales will not be forecast correctly. This could immediately 
lead to a forecast error in the planetary scales. Showalter (1984), using a spectral 
forecast verification technique, noted a case where a poor synoptic-scale forecast by 
the NOGAPS model led to subsequent degradation in the planetarv-scale forecasts of 
the model. The nonlinear linkage of synoptic and planetary scale of motion suggests 
that to forecast planetary scales more accurately, more accurate synoptic-scale forecast 


are required. 
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View haviean vex CONCLUSIONS 


The procedures used in this study are based on the normal mode analysis 
procedures developed by Errico (1984). NOGAPS analyses have been projected onto 
the normal modes and the different terms (linear,nonlinear adiabatic and diabatic) that 
affect a particular mode have been determined. The total energy for each mode and 
the energy tendencies due to the nonlinear term are determined. The effect of synoptic 
scales on the magnitude of the adiabatic nonlinear term (N,,). and the energy tendency 
due to N,, is also determined. By deleting the effect of the momentum advection terms, 
It is possible to determine the relative importance of interactions through the 
momentum or mass terms. 

The importance of the nonlinear interactions in the dynamics of planetarv waves 
has been demonstrated by computing the adiabatic nonlinear term for a two different 
data sets. The first data set is 19 days of 1ZUTC NOGAPS analvses taken everv 5 
days from 16 January 86 to 16 April 86. The second data set is a filtered version of the 
first data set. In this data set the data (u, v, T and £n p, ) are spectrally filtered by 
transforming the data to spectral space and then setting the coefficients of wave 
numbers 7-15 to zero. Reconstitution of the field then 1s a representation of the 
atmosphere without the influence of synoptic waves. The magnitude of the difference 
between the adiabatic nonlinear term computed from the original and the filtered data 
sets 1s taken as a measure of the dynamical importance of synoptic-scale interactions 
on planetary scale waves. The magnitude of this difference is on average about 20 
-30% of the total nonlinear term, although it may be as much as 60-70% of the 
magnitude of the total nonlinear term in certain vertical modes for a given dav. 

The ratio of the synoptic scale contribution (N,,,) to the adibatic nonlinear term 
of a planetary-scale mode (N,) is only one measure of the dynamical importance of 
synoptic scales. An additional measure ( the ratio of the magnitude of N,,, to the 
magnitude of the time tendency AC/At) is used to show that the ratio of the 
magnitudes of N,,, to AC/At may be large even though the ratio of the magnitude of 
Nps to the magnitude of N,, is small. This is generally true of the first three vertical 
modes, although it was especially true for the barotropic mode. For this mode, the 


linear and nonlinear terms tend to balance such that the time tendency of this mode is 
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smaller than either of the two balancing terms. Even though small when compared to 
Ny the contribution of synoptic scales may be important in determining the balance or 
nonbalance of the total nonlinear term with the linear term. In this case, Bec nmught be 
very important in determining the time tendency of the barotropic mode. 

Although the comparisons with the magnitude of N,, to the magnitude N, or to 
AC/AzZ show that N,,, could play an important role in the dynamics of planetary scale 
modes, it does not show how planetary scale modes are affected by synoptic scale 
modes. To examine the mechanism through which the synoptic scales affect planetary 
scales, the energy equation for a given mode is derived and the energy tendencies 
generated by N,, and N,, are examined. On the average, synoptic scales tend to give 
energy up to planetary scales for the barotropic 2=1 and baroclinic £=4 vertical 
modes, while energy flows from planetarv-scale waves to the £=2.3 modes. The 
positive contribution of svnoptic scales to the energy of the planetary scale 2= 1 mode 
can be a large percentage of the energy tendency generated by N,, which indicates that 
the synoptic scales play an important role in determining the time evolution of these 
modes. 

To determine how energy is being transferred, the energy tendencies generated 
by N,, with the momentum advection term deleted from the equations are calculated. 
The deletion of this term eliminates most of the transfer of energy through the 
momentum term. By examining the difference between the energy tendency generated 
by N,. With all the terms included and N,, with no momentum advections, it is 
possible to determine how much of the energy tendency due to interaction with 
synoptic scales is due to the momentum advections. It is found that interactions 
through this term are, not unsurprisingly, responsible for almost all of the energy 
transfer from synoptic scales to the planetary scale barotropic modes. Interactions 
through this term are also important for the £= 2,3 modes, but do not seem to be quite 
as important as the interaction through the mass fields (temperature advections). The 
relative importance of interaction through the mass field increases as vertical mode 
number increases. 

The results of this study confirm in part those of Gall et al. (1979). That is, 
synoptic scales do in fact play an important role in the dynamics of planetary scale 
waves and they can act to increase the energy of these waves. However, the results of 
this study also appear to conflict with the results of Gall et al. (1979). Gall proposed 


that the forcing of planetary scale waves by synoptic scale waves was mainly through 
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the planetary scale variations in the meridional heat flux convergence of synoptic 
waves, which produces a positive correlation between planetary wave upward motion 
and temperature. In this study, it is found that forcing through the momentum 
advections are also important (especially for the barotropic mode). The apparent 
conflict between this study can be explained by the differences in the data used rather 
than the mechanism. 

There is no question that the synoptic scales have an impact on the dynamics of 
planetary scales. For example, if a major cyclonic development is missed in a forecast 
it could lead to a sizable error in the planetary waves because the nonlinear effects of 
cyclones would not be properly represented. This fact means that to more accurately 
forecast planetary scale waves one must improve the forecast of svnoptic scales which 
will feed back and cause an increase in accuracy of the planetary scales. This has 
already been demonstrated to some degree because higher resolution models tend to 
give better planetary wave forecasts than lower resolution models even though the 
truncation error for the planetary waves should be negligible for both resolutions. 

It remains for future studies to determine how well numerical models represent 
these interactions and how errors in forecasting these interactions affect the planetary 
scale forecasts. Future studies might should try to examine how periods of large 
synoptic-planetary interactions are affected by changes which occur on the time-scale 
of synoptic wave. The interactions in this study were examined only every five days so 
it was not possible to see how the interactions varied with the changing synoptic 
patterns. 

Another possiblity for future study would be forecast verification using normal 
mode analysis. Using normal mode analysis, it may be possible to more fully examine 
the relationship between errors in determining the interactions between planetary and 
synoptic scales and any subsequent planetary-scale forecast error. Finally, other 
Studies have indicated that nonlinear interactions may be important in maintaining 
blocking patterns. These studies did not use normal mode analysis. Much might be 
learned about the dynamics of these blocking patterns by examining the nonlinear 


interactions using normal mode analysis. 
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APPENDIX A 


LINEARIZED HYDROSTATIC, THERMODYNAMIC AND 
CONTINUITY EQUATIONS 


The vertical structure of the NOGAPS model follows the development given by 
Arakawa and Suarez (1983). The variables are staggered in © so that all the variables 


(T, >, and V) except pressure and o are carried at the mid-point of each layer. The 


vertical structure is illustrated in Fig. A.1 


The finite difference form of the hvdrostatic equation 1s: 


bk — de+2 = Col Pere — Pr)Ona1 for k=1,3,..K —2 (A.1) 
and 

dx = os + Co(P, — Py) G, (A.2) 
Where 


al+K« al+Kx 





re a can 
a ia sk Oo (A.3) 
pPol+K pea — Pr-i 
eg emcee 
"pel +6 Bs — Prot (ed) 
and 
Gar = Angi On + Beri Fria (A.5) 
are the interpolation formulas used to produce energetically consistent equations 
where: 
Pk = O(Ps — pe) + Pe (A.6) 
2 Dk 
Ee = (A.7) 
Po 
Aire Jeb cee ape 
eo 
ey a ih (A.8) 
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Figure A.1. Vertical grid structure for NOGAPS model, 
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(A.9) 


Bee Pra — Prat 
k+1 Tee re 





PeER => Pies 
te oe = 
Ar-1= (A.10) 
ci PR — Pr-2 . 
ieee Pr Pr~3 
Be-i = om All 
Kies Koes (A.11) 
T is temperature, p is pressure and Cy is the specific heat at constant pressure. The 
geopotential at each level 1s computed using 
bx =¢,+C,(P, — Px)Ox (4.12) 
(A.13) 


¢k = Fe NS UUe OA pare = Pa}(An+i Pn)t+ 
oh Caen ee Fn 9) (ee re 


x 
w~N=k+2 


which are just integrated forms of (A.1 ) and (A.2). 
increments of 2. The above form of the hydrostatic equation can be written as: 


(A.14) 


The primed sum indicates 


Pk ia P» = ae Ceal aA 


where 0 
m=k 


G oo Col Paseo PrjAn+1 
kn — lass) ill 
Cn Pa+e = Pa)An+1 ae 
Pr 
CHUA Veh 
Pr 


m>k 


Or in matrix form as 
(A.16) 


o— ds = CT 
The finite difference form of the thermodynamic equation (Eq. 299 in Arakawa 


and Lamb:1972) in orthogonal curvilinear coordinates 1s 
—~ aan l ar 
byl nT hs + [be( FT") + bn(GT” Ik; + = Pribo( 58%, 
; (A.17) 
] Om wuAn——~ Dé ———4 
= —|(oma) OP 1% —" (ona) bem + V—(ona) bgt + NO}; 





where Cp 
Ae&An 
aay mn (A.18) 
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ee ge 
n 
Ae 
G = nu— 
m 
S=TIlo 


The overbar is a linear average in the direction of the variable indicated, 


difference taken in the direction of the subscript. 
To linearized, first subtract 




















k bo(S), 
Which gives 
] 6,(S'*. ] OTl 
k eas k 1 S89 pea. k . 
16.7;, ate Aoi Eo[S OF. — Mrmr COs By ey =(Qr)s; 
or d : - , 
TX + atta 0, _7T k+1 k= “k-1 p _F ‘s k—-1 
y+ een kl Ao. =| jean : T;| 
_ (ora), Olnt | 
Where T, is the rest-state temperature. = eo ra (Or); 


Substituting the linearized form of the continuity equation (Eqs. 
Arakawa and Lamb;1972) 


Ck4+1 = — —SF_AV- Vn don + Ona» SEO “Vn Acn Oe 
Chat = —Snai (V Vn Aon + Ohu1 81 (V > Va) don + Q_ 


éine =. —~V%(V-Valdon +Qp 


and 
T., dP, 
(o7Ta)k = ae Cy pty — 
gives 
k—2, 
TE + [ jy CET Unt (VV inden + tee Snes alY Val don 


Ac. 


aoe 


(A.24) 


166-167 in 


(A.27) 


(A.29) 


| (On-1 — 1S Sane -V)p Aon + Ok Ot (VY eae 
a . Xo 

1 OL, ae 

Pelt Gq enelY -V,)Ao, = (emi 

In matnx form 


Deo hee + T(V : Nas = (Qr)i5 


where 





(A.30) 
(on41 — 1) (Choy) Tare de 
a f(r —) k 
| |+ Men [ | Ae pr = Ack n< 
(ce+1 — 1) (G5-1). Tie, Be aan 
= a es +e =k 
es t + og [ I- Ack p,\* ar )} Aor i 
(Or41) (Gree ee ee 
——————— — | | ————— + — (71 — JI A >k 
{{ |+ ier | | Ac, P, (™ ps )} OK nr 
The formula 
dP, A Prag — Pr)| + Bree 21( Pe — Wee 
ok — Akriloe+i (Pea ~ Pe)| + Be-iloe-1( Pe = Pe-2)l , en (4.32) 
ar Ac; 
P =e k=K 
is that derived by Arakawa from the interpolation formulas (A.3) - (A.9). Note that 
the continuity equation may also be written in matrix form by defining 
Ao, 
Age (A.33) 
II = 
Ao 
so that 
6: Inc =-1175+Q, (Aga 
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APPENDIX B 
NORMAL MODES OF THE QUASI-GEOSTROPHIC MODEL 


It is convenient to use the following vertical coordinate 


Z = —In(p/po), (B.1) 
where Po is a fixed standard sea-level pressure. The vertical coordinate Z is related to 
the actual height z, and the geopotential bv the equation of state and the hydrostatic 
equation: 


az ad 
ee PO PIS, az (B.2) 


This vertical coordinate is closely related to the more familiar pressure coordinate 
system and the vertical derivatives of the two svstem are related as follows: 


g 140 
Op - paz’ (B.3) 
while the other partial derivatives are the same in both systems. The Z-velocity Z 1s 


related to ® through 


. Ww 
Z=-. i 
D (B.4) 


The basic equations of this model are: 


Op y 
—_ = RT (B.35) 
OZ 
OV - OV 
ae ; —— V =O . 
a +V rea Oe x (B.6) 
OT .OT (B.7) 
ey +V-VIP+ 227+ ZkT =Q 

AZ Cz 


We assume boundary conditions of Z= 0 at 0 and Z= Zy where Zr is Z at the top of 


the atmosphere. An expression for Z may be obtained by rewriting B.8 as 
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a ; 
WEN) Te AN 4) (B.9) 
if we multiplly (B.9) by et integrate from Z to Z7 and solve for Z we obtain 


ar 
t=¢ | e 7V-VdzZ. 
Zz 


(B.10) 
Assuming Boussinequess conditions 1.e. e°* = constant then (B.10) can be written 
° ir 
z= | V-VdZ. nil 


Next we linearize (B.5)-(B.8) about a hydrostatic and adiabatic mean state at rest by 
letting 

V=V'(z,y,Z,t); ¢=¢(Z)+¢'(z,y, Z,¢). (B.12) 
By combining the linearized versions of (B.10), (B.5). and (B.7) we obtain a single 


equation for mass which contains all of the vertical derivatives: 








oO"? Lr ri V-VdZ 
Zz ; 
OteZ x (B.13) 
The term I°(Z) is the mean state static stability and is given by 
@ foe. = - Hegaa leon 
(2) = 579g +") = ao + Haz) (B.14) 
where 
H = RT 9. (B.15) 


A vertical structure equation can now be determined from B.13 by first dividing 
by ['(Z) and then taking ¢ :éZ of (B.13) so that it can be written as 


O70 , do 
Sala fap) ee ae (B.16) 


Next We can use the technique of separation of variables to determine the vertical 


structure equation. This is done by letting 


¢ = ®, (z,y,t)o,(Z) (B.17) 


VES Vales, ton 2) (B.18) 


and substituting these expressions into (B.16). The vertical structure equation obtained 


by the above procedure is. 
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QO 10, 








a \ = —A, n° 

YANK OA ) (B.19) 
The associated boundary conditions for this equation are 

Abn _ 

Lt i (B.20) 


The vertical modes of this model are just the eigenfunction of the boundary value 
problem given by (B.19) and (B.20). Assuming I(Z) is constant then the vertical 
modes for this system are 
Z 
on = A, cos —— n=0— 
Zr i (B.21) 


These function can be used to transform the basic equation into equations for each 





vertical mode. 
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